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TECHNICAL OBJECTIVES 


1. Optimal Utilization of Laser and VLBI Observations for Reference 
Frames for Geodynamics (Grant NSG $265} 

2. Optimal Utilization of Satellite-Borne Laser Ranging System 
(Grant NSG 5265) Completed. 

3. Geodetic Utilization of NAVSTAR Geodetic Positioning System 
(Grant NSG 5265) Completed. 

4. Utilization of Range Difference Observations in Geodynamics 
(Contract NAS 5-25888) 


2. ACTIVITIES 


2.1 Investigations Related to the Establishment and Maintenance of 
a Conventional Terrestrial Reference System 


w.ll Introduction 

There seems to be general agreement that the new CTS frame 
conceptually be defined similarly to the CIO-BIH system, i.e., 
it should be attached to observatories located on the surface 
of the earth. The main difference in concept is that these 
can no longer be assumed motionless with respect to each other. 
Also they must be equipped with advanced goedetic instrumenta- 
tion like VLBI or lasers, which are no longer referenced to 
the local plumblines. 

[teller, 1981a] 

In other words, with the advent of new techniques (hence, the substan- 
tial increase of accuracy), the earth can no longer be thought of as 
rigid. This is particularly apparent considering the stated require- 
ment of one-day earth rotation and two-day polar motion averages to 
5 cm accuracy [Mueller, 1981b]. 

The following sections summarize the progress of investigations 
on the definition and maintenance of a CTS by a global network of 
VLBI and laser stations which can be thought of as defining the 
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vertices of a polyhedron [Van Gelder, 1978]. The function of the CIS 
can be divided Into two aspects. The first Is to monitor the exter- 
nal (or global) motions of the polyhedron with respect to a Conven- 
tional Inertial (coordinate) System (CIS), I.e., those motions that 
are common to all stations such as precession, nutation, earth orien- 
tation (earth rotation and poUr motion) but primarily the latter. 

The second function Is to monitor .he Internal motions (or deforma- 
tions) of the polyhedron, I.e., those motions that are not common 
to all stations such as plate tectonic motion, crustal deformation, 
tidal motions and ocean loading. Thus, the deformation of the poly- 
hedron will be defined as having no common rotations or translations 
[Mueller, 1981; Bender and Goad, 1979; Bender, 1981; Gulnot, 1981]. 
These two functions allow a convenient and practical division of the 
CTS Into an earth orientation monitoring service and an earth defor- 
mation monitoring service, although the two are Integrally related. 

2.12 Definition of the CTS 

The frame of the CTS Is defined by the adopted coordinates of 
a global network of stations at a fundamental epoch of time, t^, which 
constitutes the vertices of a fundamental polyhedron. Denote these 
fundamental coordinates by They will be estimated from an obser- 
vational campaign dedicated to this purpose In which a combination of 
several modern three-dimensional geodetic systems will participate. 

The estimation of X* by a coordinate transformation combined with a 

‘•0 

free adjustment approach Is described In [Bock, in preparation]. 

The fundamental polyhedron establishes a geometric description 
of the earth's crust at the fundamental epoch, and therefore serves 
as a reference frame for the monitoring of earth dynamics. In addi- 
tion, the CTS will include all necessary model parameters so that 
the system will be well-defined and unambiguous. These include pre- 
cession (P) and nutation (N) models to be used In connecting the 
CIT and CTS according to 
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[CTS] * SNP [CIS] 


where S will be estimated by the earth orientation monitoring service 
relative to the Initial orientation of the fundamental polyhedron at 
tQ. In addition, tidal models, plate tectonic models, etc. will be 
Incorporated In maintaining the system as well as other fundamental 
parameters such as the speed of light. 

The fundamental polyhedron Is aligned with the CI0>BIH system 
at tQ as described In [Bock, In preparation]. However, since the 
earth can no longer be assumed to be rigid, the new CTS must be 
referred to a set of axes that are fixed In some average sense In 
the deformable earth. A modified Tisserand axes definition will 
serve this purpose. It can be shown that the algorithms used to 
maintain the CTS are equivalent to this definition and therefore 
consistent with earth rotation theory. 


2.121 Tisserand axes. 


For the deformable earth, no body-fixed reference axes exist. 

In earth rotation theory :?3ually Tisserand 's axes are used. 

From the definition of the Tisserand axes, the relative angular 
momentum should be zero [Moritz, 1980] 

According to [Jeffreys, 1970] 


/// ^ ~ 

where jj is the velocity vector. 
From (1) and (2) we obtain 


/// 
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dm 

HI 
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- “x-^) 

dm 

III 

(u,-y 
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dm 

HI 

u, dm 

= 0 


III 

Uy dm 
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III 

u dm 

» 0 
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( 3 ) 



The Integrals are on the whole earth body*, thus Munk and MacDonald 
call these axes “Tisserand's mean axes of body ,*' Very roughly we can think 
that these axes are "fixed" In some average sense in the deformable earth. 
They are ideal axes and can hardly be realized in practice, since within 
the earth the velocity u is impossible to measure directly by geodetic 
means. 


2.122 Approximate realization of the Tisserand axes (ATA). 


First, we alter the volume integrals into surface integrals. This 
approximation is dynamically drastic. Because of this approximation any 
motion of the crust as a whole relative to mantle could not be detected. 

But this approximation is practically necessary for geodetic purposes. We 
call these modified axes the "Tisserand mean axes of crust ." Or roughly, 
we think these axes will be "fixed" in some average sense in the deformable 
earth crust. 

Second, the velocity vector is changed into a movement (deformation) 
vector. For an infinitesimal time interval, this change is rigorous. For 
a finite time interval, it is an approximation. 

Third, we use summation instead of integrals. The points to be sumned 
should be infinite and cover all the earth's crust. 

Last, for practical reasons, only finite points are available, and they 
are only on sol id earth surface, not on ocean. 

After these approximations, the requirements for Tisserand axes become 
n 


.E, (Ay*z - A2*y). M. 
n 

E (az*x - Ax*z). M. 
i»l ’ ^ 

n 

.E (Ax-y - Ayx). M. 


E AX. M. = 0 
i*l ^ 1 


E Ay.1 M. =0 
i«l ^ 


ih “i '*1 * ° 


a 0 
» 0 
* 0 


( 4 ) 
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These equations are used for practically defining and maintaining the new 
CTS. To make such deflneJ CT$ nieanlngful , a large number of well distributed 
stations are needed. If the stations are not well distributed and some sta- 
tions are located In areas less stable than others, the weight matrix M can be 
used. Otherwise, let H « I. 


2.123 Other proposed constraints for the maintenance of the CTS. 


Some authors proposed that the constraints for the maintenance of the 
CTS should be that there Is no conmon rotation (and translation) In the 
deformation (see [Mueller, 1981; Bender, 1979, 1981; Gulnot, 1981]). Now 
we show that the ATA also fulfills the above requirement of no common rota- 
tion or translation. 

Let be the deformation of station j at time t. 

Axi « xi - ^ (5) 

— t — t -0 


1 1 ' 

Imagine that at time t, station coordinates change from ^ to x;. 

i * i -0 — c 

Through a translation and a rotation Is changed into xj, so that 

In other words, by doing 


lit -2o 


J|2. 


min, where w. is the “weight." 


i — V -V j 

so, x^ will have no cotranon rotation or common translation with respect to 
)(q. Let R(^) be the rotation and 6x be the translation. Then 



4 • 

R( 6 .) j - Ax 

(6) 


“t • 

^ • R(9)*X^ - fix - 

““C —0 — 

(7) 

Form 

s » 

I Wj lixJp 

(8) 

By 

- 

35 “ 

Hx ■ « 

(9) 


we obtain the six constraints. 


In a Cartesian coordinate system, equation (9) becomes 
z Wj*AXj » z Wj aYj » 2 Wj aZj » 0 (10a) 

z Wj(aY‘Z - AZ*Y)j • iWj(aZ-X - aX-Z)^ - zwj(aX-Y - aY*X)j. - 0 (10b) 


Equations (10) are exaccly the same form as equation (4) for ATA require- 
ments. 
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For the geodetic coordinate system (H, X), the above equation 

changes Into 

£ Wj(*R cos^ sinx ax • R sln^ cosX cos^ cosx aH)j ■ 0 

£ Wj(R cos^ cosx ax > R sln^ sInX a^ co$<^ sinx aH)j ° C (11a) 

£ Wj(R cos^ a^ '** sln^ aH)j ■ 0 

£ Wj(cos^ sln^ cosx ax - sinx a^)j « 0 

£ Wj(cos^ slni^ sinx ax + cosx a(^)j > 0 (lib) 

£ Wj(cos^t ax)j « 0 

For the local horizon coordinate system where x points east, y north, 
z upward, neglecting the earth's flattening, the constraints are 

£ Wj(“Slnx ax - sii<^ cosX ay + cosa cosX az)^ * 0 

£ Wj(cosX ax - s1.<^ sinx ay * cos^ sInX az)j « 0 (12a) 

£ Wj(cos^ Ay + sln^ az)j ■ 0 

£ Wj(s1n^ cosx ax - sInX ay)^ ■ 0 

£ Wj(s1n^ sInx ax + cosx Ay)j “ 0 (12b) 

£ Wj(cos^ ax)j * 0 

Equations (12b) are the same as given by Bender and Goad [1979]. 

These constraints also require a sufficient number of well -distributed 
stations. It is hoped that there Is no cornnon rotation (or translation) with 
respect to the whole earth crust. If there are not enough stations and/or 
they are ill distributed, then these station deformations do have some cornnon 
motion. Using these constraints to eliminate the common motion In the defor- 
mations will bring the true cornnon motion Into the determined polar motion. 
That Is, we cannot separate the stations' common motion from real polar 
motion. 

When a sufficient number of well -distributed stations are available, 
the CTS defined and maintained by these contraints are ATA, and at the same 
time It fulfills the condition of “no common rotation and no common trans- 
lation." 
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2.124 Dimension of the new CTS. 

The CTS's of different observing techniques have different dimensions, 
depending on the observables and also on the coordinate system of the observ- 
ing targets (CIS). 


Techniques: 

Optical 

VLBI 

LR, GPS, Doppler 

CIS: 

2-d1mens1onal on 
unit sphere 

Same as optical 

3-d1mens1ona1 

Observables: 

Angle 

Distance (time) 

Distance (time) 

CTS: 

2-d1mens1onal on 
unit sphere 

Quasi -3-d1mens1onal 

3-d1mens1onal 


For VLBI the observable Is distance so that the CTS built by VLBI Is three- 
dimensional, but the source position is two-dimensional (on a unit sphere). 
Therefore, the VLBI CTS Is not fully three-dimensional. One can only deter- 
mine distance and shape, but not the center of the coordinate system. We 
call It quasi -three-dimensional . The new CTS seems to be the combination of 
a few new techniques. Therefore the new uniform CTS Is three-dimensional. 

It seems insufficient to define a three-dimensional new CTS only by longi- 
tude origin and pole position. 


2.125 Sensitivity of the rotation with the height change. 

In a Cartesian coordinate system, a small rotation ^ and a translation 
^wlll change the coordinates In such a way 

■ 03 Y^ - 62 !^ + 6 X 

* 6 iZ^ - fljX^ ♦ 6 Y (13) 

• 02X^ - flxY^ + 6Z 


The 

Ah^ 


corresponding height change Is 


aZ^-Z^ + aY^-Y^ + aX^.X^ 


(14) 


If only a small rotation exists, that Is, £ » 0, then we will find 
Ah^ a 0 


(15) 
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It means that the rotation change only causes horizontal coordinate changes. 
The height Is not sensitive to rotation change. The above result Is exactly 
as stated by [Bender and Goad, 1979]. 

But the height being Insensitive to rotation does not necessarily 
mean that rotation Is also Insensitive to height change. On the contrary, 
we now prove that rotation Is sensitive to height changes (In the practical 
station distribution). 

Now the question Is to solve for e and ± from the height change by 
equation (13). In this case a least squares method can be used. The normal 
equation Is 

N X B * U (16) 


where 



U'+Y*] -CXY] 

-CXZ] 

0 

[Z] 

-CYf 


[X^+Z^] 

-CYZ] 

-CZ] 

0 

[X] 



[X*+Y^3 

[Y] 

-[X] 

0 

N « 



n 

0 

0 


symmetric 



n 

0 


m 




n. 



ei* 


Iaz-y - AY-Z]' 


02 


[aX-Z - AZ-X] 

8 - 

01 

U * 

[aY‘X - AX-Y] 


[AX] 


«Y 


[AY] 


.5Z, 


. [AZ] 


For the practical station distribution, the nondiagonal elements In N 
generally are not equal to zero. We write N’^ as 


Now we transform the above equation from the Cartesian to the 'ocal horizon 
system, and letting > ax > 0, then 

0 

0 

0 


u 


[Ah cos^ cosx] 
[Ah cos^ sinx] 
iJ[Ah sln^] 


(18) 
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(19) 


®i ®i4 cosx] + a^g[Ah cos<j> sinx] + a^g[Ah sin(j»] 

i = 1.2,3 

Since a^g, a^g generally are not equal to zero, 0^ also is not equal 
to zero. Simulation shows that can be of the same order of magnitude 
as Ah. That means that rotation is sensitive to height change. 

In equations (10), (11) and (12) we deliberately divide the constraints 
Into two groups: (a) for translation, (b) for rotation, although for practi- 
cal reasons people seem to be more Interested in rotation than in translation. 
However, through the above derivation we see that height change will, through 
translation, influence rotation. Since the CTS is three-dimensional, the six 
constraints are a complete set, and we could not use only set (b), even if 
we are only Interested in rotation. 

2.13 Maintenance of the CTS 

2.131 General Description. 

The primary functions cf the CTS are twofold. First, to moni- 
tor the variations in earth orientation relative to the fundamental 
epoch and second to monitor earth deformation. The first function 
will require a dedicated network of observations on a continuous basis 
considering the requirements of one-day earth rotation and two-day 
polar motion averages to 5 cm accuracy. The second function does hot 
require such an intense observational schedule since it Is antici- 
pated that the time variation of the polyhedron coordinates due to 
deformations will be in the submeter range per year considering 
plate tectonic theory [Minster et al., 1974, 1978]. In addition, 
solely due to economic and other practical considerations, we could 
not expect all the CTS stations to observe continuously. Further- 
more, the first function requires dedicated fixed observations while 
the second function may also involve incorporating mobile stations. 

Vie propose the following setup for the CTS as Is depicted in 
Fig. 1. 
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Fig. 1 Schematic CTS operations. 
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The polyhedron Is composed of a set of stations distributed (In a 
manner described In [Bock, In preparation]) globally. In the figure, 
we see an observation schedule divided Into two total Intervals, 
A(tx-to) and B(t2-ti). Level I, the Earth Orientation Monitoring 
Service, Is composed of a dedicated subset of the GTS. These sta- 
tions monitor earth orientation on a regular basis denoted by the 
subintervals within the total Intervals. These results are made 
available to the user community as Is done by the earth orientation 
services today. Level II, the Earth Deformation Monitoring Service, 

Is composed of all the polyhedron stations which at the end of each 
total Interval observe together In one short campaign (denoted In the 
figure by the shaded portion). It seems evident that observations 
from a larger number of stations, distributed over several of the 
larger tectonic plates will provide a better representation of the 
deformation of the earth than would the observations from the limited 
(due to economic and practical considerations) number of stations 
of Level I. From the Level II solution which will basically have 
the same parametrizatlon as Level I, the baseline chord lengths will 
be extracted along with their covariance matrix. These will be used 
as Input for a generalized free adjustment solution (to be described 
in [Bock, In preparation] to estimate corrections to the polyhedron 
coordinates due to deformations. For example, as shown In the figure, 
at the end of Interval A, the Level II solution will yield 
which will be added to the fundamental coordinates, to be used as 
Input In Interval B. In this way, we prevent the Internal deforma- 
tions of the polyhedron from contaminating the global rotation param- 
eters and we keep the GTS referred to a set of axes that are fixed 
In an average sense In the deformable earth. The procedure Is repeated 
for each total Interval in the same manner as depicted In the figure. 

It should be mentioned that the polyhedron coordinates may also be 
corrected, as shown In the figure, by Vj^(t-to) where v^^ (the time 
variation of the coordinates) is derived from a geophysical model 
adopted by the GTS, If an accurate enough model Is available. This 
aspect Is discussed In [Bock, In preparation]. 
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2.132 Level I - Earth Orientation Monitoring Service 

The following discussion Is restricted to VLB! observations 
since it Is felt that the Earth Orientation Monitoring Service will 
be composed primarily of dedicated VLBI observatories. This measure- 
ment system has distinct advantages over the other techniques. VLBI 
measurements are Independent of the gravity field of the earth making 
them essentially geometric In nature. Furthermore, the system has 
virtually all weather capability. Most Important, VLBI has greater 
sensitivity to earth orientation parameters than do satellite laser 
observations particularly In UTl-UTC variations [Van Gelder, 1978]. 

The VLBI mathematical models for delay and delay-rate observa- 
tions are described In detail In [Robertson, 1975; Ma, 1978]. For 
the purpose of this discussion we will follow the notation of [Bock, 
1980]. 

The choice of parameters is particularly important in the earth 
orientation solution (Level I) since all earth orientation parameters 
need to be related unambiguously to the initial orientation of the 
fundamental polyhedron given by which has been made to coincide 
with the CIO-BIH system at the fundamental epoch, Lq. As shown in 
[Bock, 1980], the earth orientation parameters Uo* Hq* <q) and the 
baseline coordinate differences are inseparable since the initial 
orientation of the terrestrial frame is not sensed by the observables. 
In practice, this dependency is broken by not parametrizing earth 
orientation over, say, the first day of observations of a particular 
observation campaign. Subsequently, three earth orientation param- 
eters are estimated over so-called earth orientation steps, each step 
spanning a certain period of time [Dermanis, 1977]. Thus, the earth 
orientation parameters average values over 

each step relative to (Sq, Hq, Kq) '^’^Ich is input into the adjustment 
and used over the first step. However, as shown in [Bock, 1980] any 
errors in the initial orientation assumed for the first step biases 
the estimates of coordinate differences (or coordinates if in addi- 
tion one station is fixed as is usually done in practice) making 
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this parametrlzatlon particularly unsuitable for the maintenance of 
a CTS. It is, furthermore, difficult to relate the orientation ori- 
gin of different data sets. In addition, the fixing of the first 
step Is done In a subjective manner and the resulting parameter esti- 
mates are Influenced by the length of the step and distribution of 
observations over that step. 

The only geodetic parameters estimable from a polyhedron of 
VLBI stations are the size and shape of the polyhedron, i.e., the 
baseline chord lengths and the orientation of the polyhedron with 
respect to the fundamental polyhedron given by three rotation 
angles. Thus, by parametrizing baseline chord lengths and rotation 
angles we can avoid the bias introduced by fixing a step since In 
this case It is no longer necessary to fix one. After all, baseline 
lengths are Invariant with respect to coordinate system definition. 
By Introducing as Input into the adjustment X. (corrected for poly- 
hedron deformation), the three rotation angles (c* n, k) per step 
always refer to the fundamental orientation of the polyhedron, unam- 
biguously. 

The partial derivatives of the design matrix for the VLBI ad- 
justment are given In [Bock, 1980]. The contribution to the design 
matrix for baseline chord length, B.., is given by 

" J 

where A denotes the partial derivative of the observables (delay) 
with respect to the subscripted parameter (see [Bock, 1980; Bock, 
in preparation]). 

It should be mentioned that the baseline chord lengths propa- 
gated from the estimation of coordinates (or coordinate differences) 
are unaffected by errors in 5 ^, tc^ fixed over the first step. 
However, this type of adjustment is less rigorous because of an 
arbitrary fixing of the first step to provide the necessary minimal 
constraints. Our proposed solution involves only the strictly esti- 
mable quantities, I.e., there is no need for any constraints except 
for at least one source right ascencion to fix the origin of the 
inertial frame. However, any error in this value, dao, will not 
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propagate Into baseline lengths or polar nwtion but only Into k. 
the earth rotation parameter. By being consistent with the origin 
choice (ao), 1.e., adopting a fundamental source catalogue (CIS), 
this error Is essentially eliminated. As far as parametrizing radio 
-source coordinates (a, 5), It is probably best to hold these fixed 
according to the fundamental coordinates of a source catalogue, and 
updated to the time-of-observation epoch. 

2.133 Level II - Earth Deformation Monitoring Service 

In the analysis of earth deformation, all CTS stations observe 
according to a specified schedule, to be described in [Bock, In prep- 
aration], the point being that the observational period is short and 
infrequent. In contrast to Level I, the emphasis is on the estima- 
tion of baseline lengths (rather than earth orientation) as a key to 
the analysis of the deformations of the polyhedron with time. Thus, 
any of the systems that provide accurate baseline lengths can partic- 
ipate. 

The estimation of earth deformation proceeds in two steps. 

First, the observations of the Level II campaign are analyzed as in 
Level I except that the Input prarameters for the station coordi- 
nates are held at , possibly corrected for the motions inferred 
^0 

by an adopted geophysical model as described in [Bock, in preparation]. 
In this way, the connection to the fundamental polyhedron is maintained. 
The parameters that are estimated in the first step are the same as 
in Level I— baseline chord lengths, earth orientation angles (in 
order to remove all cocmion rotations) and model parameters. 

For the second step, the estimated baseline lengths as well as 
their estimated covariances are extracted from the first step adjust- 
ment. From this information, the size and shape of the deformed 
polyhedron is completely defined, assumed constant over the short 
observational campaign. However, the absolute location of the poly- 
hedron, i.e., its new coordinates are undetermined by the change of 
its sides, but this is exactly what we seek. The problem is singu- 
lar, being rank deficient by 6, the familiar origin and orientation 
defects. In order to overcome this singularity we will employ a 
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free adjustment approach generalized to Include an a priori geophys- 
ical model for the time variations of the polyhedron vertices. It 
will be shown that this method insures that the CTS axes are fixed 
in the sense of the modified Tisserand axes definition, in the de- 
formable earth. Furthermore, it will be shown that the estimation 
method yields parameters with optimal statistical and geometric 
properties. This approach has been applied to crustal deformation 
analyses (see [Brunner, 1979; Brunner et al., 1980]) and has been 
suggested by [Dermanis, 1980; Bender and Goad, 1979; Bender, 1980] 
for application to terrestrial reference frames although in a less 
general and detailed manner than will be presented in [Bock, in 
preparation]. 

2.134 Some comnents about the practical realization of the new 
CTS. 

For some conventional reason, we want the new CTS to be the con- 
tinuation of the BIN system. Let tQ be the epoch for alignment. This 
means that at epoch to the z-axis of the new CTS will coincide (with- 
in the scope of error) with the BIH-CIO pole. Therefore, the new 
CTS generally cannot be Tisserand axes at time tQi even in the approx- 
imation point of view. By applying the six constraints we keep the 
relative position between the CTS and the Tisserand axes at any time 
t the same (approximately) as at time to> so that the polar motion 
(earth rotation) in CTS is the same as that in the Tisserand axes. 

As stated before, we must have a large number of well -distrib- 
uted stations to make these constraints meaningful. By deformation 
we generally mean the secular or long-period deformation. For econ- 
omical reasons, we could not ask all the CTS stations to observe con- 
tinuously through the year; only once in a while is necessary and 
possible. On the other hand, the earth rotation service stations 
need continuous observation. Therefore, we can select a number of 
CTS stations to be the service stations. Since these service sta- 
tions' deformations are determined by the process of CTS maintenance, 
their distribution only should be optimal for earth station parametar 
determination. 
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2.135 Ongoing Investigations. 

The algorithms for both Level I and Level II are being tested 
by simulations as well as with real data. An error analysis will be 
performed for earth orientation and deformation parameter estimation. 
In addition, optimal design questions such as CIS station distribu- 
tion, observation regularity and scheduling will be investigated. 
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2.14 Comparison of Polar Motion Data from the 1980 Project MERIT 
Short Campaign 

(presented at lAU Colloquium No. 63 on "High Precision Earth 
Rotation and Earth-Moon Dynamics: Lunar Distances and Related 
Observations," May 22-27, 1981, Grasse, France) 


1. INTRODUCTION 

The short campaign of MERIT (a program of international collaboration 
to Monitor E^arth Rotation and ^ntercompare the Techniques of observation 
and analysis) was held during the three-month period August to October, 
1980. The participation was quite extensive and varied in techniques: 82 
instruments from 22 countries provided classical astrometric data; the U.S. 
Defense Mapping Agency and GR6S (France) through the MEDOC network supplied 
Doppler satellite observations; 900 passes of Lageos and 780 passes of 
Starlette were observed from as many as 27 laser ranging stations and ana- 
lyzed by SAO (Smithsonian Astrophysical Observatory, Cambridge, Mass.), 

CNES (Centre National d‘Etudes Spatiales. Toulouse, France), UTX (Univer- 
sity of Texas at Austin), and 6SFC (NASA Goddard Space Flight Center, 
Greenbelt, Md.). Comparison of the above data set is the subject of this 
paper. Additional observations were made by connected-element radio inter- 
ferometry at Green Bank, W.Va., and at Cambridge, U.K.; very long baseline 
radio interferometry at about ten stations around the world; lunar laser 
ranging at the McDonald Observatory, Texas. Data from these techniques 
is not included in this comparison because they either did not provide con- 
tinuous data throughout the campaign, or the data was not readily avail- 
able through the Coordinating Center at the 81H in Paris. 
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2. PREDICTION OF POLAR MOTION 

Based on the analysis of polar motion behavior, the possibility of 
predicting polar motion for a long time Interval (1-2 years In advance) 
with sufficient accuracy has been found. The best estimated Chandler 
period Is taken as constant, and six years of data are used to estimate 
the amplitudes, phases and elllpticity of the Chandler and annual motions. 
These estimated parameters are then used tu predict the next year's (or 
next two years') polar motion. In making the prediction, possible linear 
trends are also taken into consideration. 

The data used for prediction were those of the 6IH, IPMS and DMA, 
from 1968 to 1980 (DMA from 1972 to 1980). Polar motions have been pre- 
dicted for one year In advance and compared to smoothed observed ones: 
the mean rms of the differences (predicted minus observed) Is about 01*02. 
Differences between relative polar motions are much small er: for a time 
Interval of 20-30 days, the rms difference Is about 0'.'01 (30 cm) through 
the whole year. Compared with the best available VLBI results (from 1977 
to 1980), the rms difference Is only 0U013; and the rms difference of the 
relative polar motions (with time Interval less than or equal to two 
months) Is 0'.'008, both values being remarkably small. 

The predicted polar motion can be used for geodetic purposes. It 
seems that the accuracy of the prediction Is high enough for any practical 
purpose that requires real time polar motion up to an accuracy of, say, 

50 cm. This would Include rather sophisticated applications such as con- 
trol of space probes. Details may be found In "Prediction of Earth Rota- 
tion and Polar Motion," Dept, of Geodetic Science Rep., Ohio State Univ. 
(In press). 


3. MERIT DATA INTERCOMPARISON BASED ON RAW DATA 

The BIH has already analyzed the MERIT data using BIH Circular 0 as 
a reference. Since new techniques are expected to obtain better accu- 
racies than now available from the BIH, It was thought that a mutual com- 
parison may be more rational than using Circular D as a common reference. 
Both raw data and smoothed data were used. 

3.1 Absolute Comparison 

The standard errors of the polar coordinate differences of each two 
Analyzing Centers were computed In all combinations with the following 
results: 
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DMA ~ MEDOC 

CNES(Lageos) - CNES(Star1ette} 

SAO * CNES(Starlette) 

Astrometry - SAO 

Astrometry - DMA 

Astrometry - UTX 

DMA * SAO 

UTX - SAO 

UTX - DMA 

SAO • CNES(Lageos) 

UTX - CNES(Lageos} 

SAO - Circular 0 
UTX - Circular 0 



®Ay 

mean 

0?040 

0?044 

0?042 

0.031 

0.053 

0.042 

0.040 

0.036 

0.038 

0.022 

0.025 

0.024 

0.022 

0.025 

0.023 

0.019 

0.026 

0.023 

0.010 

0.011 

0.011 

0.014 

0.012 

0.013 

0.014 

0.010 

0.012 

0.009 

0.014 

0.012 

0.010 

0.017 

0.014 

0.009 

0.009 

0.009 

0.011 

0.013 

0.012 


From the above values it is suggested that the standard errors are 
approximately 0'.'04 for MEDOC and CNES(Starlette); 0*.'02 for classical 
astrometry; 0'.'008 - OVOl for SAO, DMA, UTX and CNES(Lageos). 

3.2 Relative (Variation of Polar Motion} Comparison 

Relative polar motion comparison is used to detect possible system^ 
atic errors (other than a constant due to the difference of origins). Sys 
tematic errors will most likely be smaller in relative polar motion. 

The standard deviations of polar motion variations are given below. 



®6AX 

®6Ay 

mean 

DMA - MEDOC 

0'.'039 

0^039 

0'.’039 

UTX - SAO 

0.010 

0.010 

0.010 

UTX - DMA 

0.014 

0.008 

0.011 

DMA - SAO 

0.010 

0.009 

0.010 

Astrometry - SAO 

0.017 

0.022 

0.020 

Astrometry - DMA 

0.024 

0.022 

0.023 

Astrometry - UTX 

0.017 

0.024 

0.021 

SAO - Circular D 

0.005 

0.004 

0.005 

UTX - Circular D 

0.011 

0.010 

0.010 


Since every o^x is more or less larger than og^x* appears that there 
may be systematic errors among the Analyzing Centers. 


4. MERIT DATA INTERCOMPARISON AFTER SMOOTHING 

The above comparisons gave reason to expect systematic differences 
between the results of the earth rotation parameters as obtained from dif- 
ferent techniques. A sR»othing process was attempted to find the proper- 
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ti es of t hese differences In ter ys of th e amplitudes of the Chandler 
+ k\) and annual motions (AJ ♦ ki) and the centers of the polhodes 
(ki and k*), using the following well known circular model: 

X * kx + ka cos A + ks sin A + k% cos C + k« sin C 
y * k« - ka sin A k| cos A - k% sin C ks cos C 

where A ■ 2ir (MJO - 42413 )/365 (annual frequency) 

C > 2ir (MJO - 42413)/435 (Chandler frequency) 

These coefficients are listed in Table 1 together with their predicted 
values obtained as explained In Section 2. The values of coefficients ki 
and ks were plotted on a graph with their standard deviations (see figure). 
They correspond to the x and y coordinates of the centers of the circles 
depicting the pole movement. This figure shows that systematic differ* 
ences may exist In the pole origin. The tabulated amplitudes of the an- 
nual and Chandler motions Indicate the same thing. The agreement between 
the predicted and SAG values Is truly remarkable. 

The relatively large standard deviations and the shortness of the 
data span naturally cast a shadow of doubt on the validity of these coef- 
ficients. For this reason another adjustment was performed In which the 
amplitudes of the Chandler and annual motions (coefficients kz - ks) were 
constrained to their predicted values and only the coordinates of the pol- 
hode centers (kx and k$) were computed. The results are In the last two 
columns of Table 1 and In view of the preliminary data set the agreement 
Is remarkable. There seems to be very little evidence of systematic dif- 
ferences In the pole origin. 


5. CONCLUSION 

The data available for the MERIT Short Campaign as expected was not 
really sufficient to enable us co arrive at conclusive evidence regarding 
the systematic differences and their estimation. The noise level was also 
quite high; thus the Importance of the Intercomparison of data to be made 
available through the MERIT Main Campaign cannot be overemphasized. 

All results In this paper were based on the data distributed during 
the campaign by the Coordinating Center at the BIN except for the GSFC 
data. 
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2.2 Utilization of Bange-Difference Obser?ations 
in Geo4ynanics 


2.2i utilization of Sinai taneooslaaeos Ra noe-Diffetencea 
in G^odTnanics 


introdoction 

h snanary of the research coapleted under this project ^ was 
presented in the 4th Lageos VorluLng Group fleeting, held at 
BASh-GSfC ('^cpteiber 1t2,1981). The presentation is susuarized 
in Appendii ' « In the following sections we elaborate on sone 
theoretical <^nescions that we feel to be of najor iwpoctance in 
understanding the need for developing the sinnltaneons 
range-differencing technigne. A nunber of prelininary results 
obtained fros a very snail data set are also discussed. 


2.211 Bstinability of the polar notlcMi step function. 

It has been well established by now (e.g., [Yan Gelder, 1978 ]) 
that the coordinates of the pole are inseparable fros the 
observing static positions; a change in one can be attributed to 
either of the two. He therefore seek for a way to circus vent 
this deficiency and uncouple the two paraseters. In [ibid.], the 
idea of detersining variations in the notion of the pole has been 
proposed, instead of absolute polar coordinates. Even in this 
case however, the variations are only "conditionally estinable," 
the reason being that an origin, with respect to which the 
variations are referenced, sust still be enforced and obviously 
the bias on this initial position affects the estisated 
variations. Fortunately, translational biases do not contasinate 
the solution, but a rotation of the coordinate systes (coning 
fros errors in OTI-OTC values) does affect the pole variations. 
This error is expected to be insignificant though, considering 
that even with present capabilities the errors in UT1 are no 
greater than 1 ns . He assuae that this bias can be tolerated and 
we proceed to establish a sethod by which the initial position 
inforsation will be introduced in the estisation process. 


2.212 Discussion of the estisation process. 

There are veirious ways of introducing the definition for the 
systes of coordinates of the pole, each having its own 
intricacies, while still sain tain ing the continuity with the 
present BIH-1979 systes. He could rotate the station positions 
for instance into the sodified terrestrial systes which has its 
Z-azis coincident with the sean position of the instantaneous 
pole over the first step or its position at the beginning of ohr 
caspaign. He assuae the coordinates of the stations fixed to 
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these valaes (except for systematic effects* o*9«# tides*etc. for 
which we correct a priori) and we therefore only solve for the 
coordinates of the pole which are nothing else but the variations 
with respect to the adopted initial pole position. This nethod 
is sinple in concept* but has the distinct and undesirable 
disadvantage that any biases fron the station coordinates will 
propagate into the pole position directly. 

On the other hand* we could solve for the positions of the 
pole and the ground stations simultaneously (except for an origin 
of longitude definition)* provided that we assume the position of 
the pole over the initial step is perfectly known. In this 
sense* the variations of the pole at each subsequent step will 
refer to this initially adopted position. The station 
coordinates which are being adjusted sinultaneously* must be 
primarily determined from the observations on the initial step. 
This ofcourse indicates that the duration of the initial step as 
well as the amount and dxstribution of the data over that period* 
will be crucial in defining the initial orientation of the ground 
network. It goes without saying that long gaps in the data 
record cannot be tolerated in this scenario* since the quality of 
the solution strongly relies on a geometrically sound 
satellite -observer configuration. 

It is obvious by now that since we are using 
quasi-simultaneous range-differences which reduce by a large 
percentage the anount of usable data and that* in addition to 
this* optimal configuration requirements further reduce this 
amount* we have to have very high quality data or the accuracy of 
our estimates will be very poor. It is anticipated therefore* 
that if this technique is adopted for continuous monitoring of 
the notion of the pole* the observing schedule of the 
participating stations should be coordinated in such a way that 
all possible Lageos passes are coobserved when possible and that 
redundant baseline pairs are included in the network to 
coapensate for data loss froa poor weather* station breakdowns 
and satellite observability problems. 

The two approaches ontlined above are quite sisilar in the 
way that the origin and orientation of the polar coordinate 
system is defined. In the first case the bias in the definition 
will come froe the adopted initial position of the pole pins any 
unaccounted for station notions since we are not solving for 
their positions, in the second case* however* the coordinates of 
the stations are adjusted to fit the data as well as the 
initially adopted pole position. It is therefore obvious that 
any unfavorable data distribution or poor data gnality during the 
initial step* will directly affect the estimated pole variations, 
aoreover* the station coordinate system definition changes fron 
one solution to the next* in a rather uncontrollable manner* so 
that not only the orientation of the two systens is biased 
(provided by the estimated pole variations) * but the origin 
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definition itself. 

Since none of the above procedures will produce the level of 
accuracies that we are after, we eust devise a different approach 
where the effects of the aforeaentioaed control factors will be 
either eliainated or at least tept below tolerance levels. The 
aethod proposed here is the following. 

1 set of consistent station positions for the aonitoring 
observatories is obtained initially froa a long arc (5-6 years) 
of Lageos data. The sienltaneously obtained satellite epheaeris 
defines a guasi-inertial systea (dynaeical CIS), which is linhed 
to the terrestrial systew via an adopted set of precession, 
nutation and polar notion aodels. the latter being possibly the 

systee so that continuity is preserved. Once this master set 
of station positions and the satellite epheaeris have been 
obtained, the changes of the station positions are monitored with 
respect to the evolving satellite epheaeris (for which we do not 
adjust anymore) . The discrepancies in the adjustaent are 
attributed to the variations of the position of the pole. The 
advantage of this systea over the ones previously described is 
that in this case the original relationship between the inertial 
and the terrestrial systea is being preserved within the bounds 
of observational accuracy. Purtheraore. since the terrestrial 
systea has been defined on the basis of a long data record, any 
biases coning from the polar notion sodel used in the reduction 
should average out. 


2.213 Orbital bias elimination by range-differencing. 

The force aodel that governs the notion of a satellite is not 
perfectly Xnown and so its inadegnacdes vill contaninate the pole 
positions. If though guasi-sinultaneous ranging is used, the 
orbital aodel inadequacies vill have little effect on the 
results. In fact to understand how this is achieved, let ns 
assuae that onr data set consists of tvo subsets of ranges where 
each range in the first set L . has been obtained siaultaneously 
with the corresponding range in the second set L • The 
range-difference data set is therefore : 

L » Li - Li (^) 


Based on the satellite epheaeris and station positions, a set of 
conputed range-differences is obtained. The conparison of these 
to the observed ones during the estiaation process vill. 
of course, result in adjusting the pole position to fit the data 
and therefore estinate its notion, issune now that for each 
position of the satellite, the orbital aodel introduces a bias b 
in the conputation of the range. This bias is not the sane at 
each instant but this need not concern ns for the tine since the 
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ranges are dittecenced only at corresponding epo<^s« The 
conpated function values (ranges) can be separated into teo 
coniK>nent.s; one uhich is the value they would have if the bias 
were not there, • and the bias itself Cb^^b^) : 

l>i • lI + bi (2) 


and 

Lj * lJ + 62 (3) 


The equivalent for range-dif terences then is : 

- Li - Li («) 


or 

^RD * lI - lI + bz - bi (5) 


If the bias for each pair of ranges were the sane, then b «b 
and the resulting conputed differences are bias-tree. 
Onrortnnately though, this cannot be assnned since different 
Bodel deficiencies have different effects on the satellite orbit 
and then depending on the relative location of stations and 
satellite positions these biases propagate in a different Banner 
in the conputation. 


To Bininize these biases one possibility is to apply the 
differencing one wore tine. By differencing consecutive 
range-differences these biases alnost disappear since their 
behavior cannot change significantly over the short interval 
(a few seconds) between two consecutive observations. So if for 
the ith range-difference : 


4. • ’■li ■ ■ •’'ll 

and for the ith-fl : 

* 4i+l ‘ ^21+1 

then the doable difference is : 

*t)D^ 


( 6 ) 


(7) 


( 8 ) 
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or 

^DD. “ ^^2i+l “ hl+l^ ’ - (b2-bi)^] (9) 

Th« tera in the brackets will be extreeely saall and it is 
therefore hoped that the final aisclosare will be coapletely free 
of the orbital aodel biases. Care should be taken to coapate 
independent doable differences (that ist use obserwations 1'~2# 
3-4, 5-6, etc. and not 1-i, 2-3, 3-4, etc.) . If correlated 
observations are used then the correlations should be incladed in 
the error propagation process or the adjnstaent results are 
aeaningless. In this case however, the coapatational burden is 
such that it aakes the problem unsolvable in practice. 


2.214 Preliainary results from the 1979 Lageos data. 

The laser systeas operating daring this tine period were not 
particularly well distributed for providing sinultaneous tracking 
of Lageos, and even worse, the ones that had collected any 
significant aaounts of data were all on the continental QSk (the 
five laser-TLfil intercoaparison campaign stations) • It is not 
surprising, therefore, that only a few short tine intervals 
(3-5 days) existed in this data set which could be used to test 
the developed software. The November 27 through November 
30,1979, interval was finally selected as the one to be used for 
the tests. Stations which had sinultaneous tracking during this 
period are : STAI.AS at (iSFC(7063), Haystack (7091), Owens Valley 
(7114), Goldstone (7115), and Ft. Davis (7086). It is worth 
notinq that although about 21,000 ranges were collected during 
this period, only 1/7 of thea fell in the sinultaneous tracking 
category. Out of these 3000 ranges a data set of about 1000 
range-differences was coapiled using the technique outlined in 
Appendix 1, taking into account that the resulting differences 
should be uncorrelated (assuming that the original ranges are 
such) . The statical locations and the siaultaneous events on the 
six Lageos passes are shown in Fig.1 . In Tables 1 axid 2, the 
nuaber of events per station and per baseline respectively are 
given tor each pass. The totals in the bottom' line of the first 
table indicate the nuaber of independent events in each pass, 
while the colunns indicate how many tines each station 
participated in an event. The grand total is the total nuaber of 
events in the experiment. The numbers in the second table have a 
siailar aeaning, except that we are now considering station pairs 
rather than the stations independently. 
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Table 1 sianltaAeous Zvents DistrxbQtion by Station 


Station 


liSgeos 

Pass 



Total 

1 

2 

3 

4 

S 

6 

7063 

206 

0 

108 

14 

108 

52 

488 

7086 

5 

0 

0 

73 

108 

60 

246 

7091 

0 

155 

173 

113 

111 

19 

571 

7114 

201 

0 

106 

0 

0 

11 

320 

7115 

0 

155 

173 

54 

111 

0 

493 


Total 206 155 281 127 219 71 1059 


Table 2 sianltaneous Events 

Distrlbntion 

by Baseline 




Lageos 

Pass 



A 1 

Oo S^XXu8 

1 

2 

3 

4 

5 

6 

£OZm1, 

7063-7114 

201 

0 

108 

0 

0 

0 

309 

7063-7086 

5 

0 

0 

14 

108 

52 

179 

7115-7091 

0 

155 

173 

54 

111 

0 

493 

7091-7086 

0 

0 

0 

59 

0 

8 

67 

7091-7114 

0 

0 

0 

0 

0 

11 

11 

Total 

206 

155 

281 

127 

219 

71 

1059 


Ezaaining these tables, and especially the second one, it is 
obvioQs that only the first (70 63 -7114) and the third (7091-7115) 
baselines have any significant nonber of observations to produce 
a eeaningfttl estinate. In fact the two together account for nore 
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than 3/4 of the total oaaber of events. Qnfortaiiatel 7 « loolting 
at Pig. i it is evident that all of the observed passes are in 
the sane general direction. The two baselines sentioned above 
are also so nearly situated that the fact that we have data fros 
both is irrelevant; we would have the saae asount of inforsation 
if we had the sane nuaber of events fros either one only. So due 
to this poor geosetry of our network it is iaprobable that any 
biases will be elisinated by sisple differencing. In aost cases 
the range froa one end of the baseline is several tiaes longer 
coapared to the one froa the other end, and in a quite different 
direction with respect to the satellite groundtrach. 


2.215 Besnlts and conclusions. 

Tn light of the above reaarhs one would therefore expect no 
astonishing results out of this experiaent. As a natter of fact 
these tests were done to check the developed software rather than 
the underlying aethod. Certain conclusions though can be still 
safely drawn, so that we consider then in future ezperiaents. As 
far as polar notion is concerned, a well balanced set of passes 
is very iaportant. Although a pair of orthogonal baselines can 
provide estiaates for both of the conponents in theory, a 
worldwide set of such evenly distributed pairs is what we should 
be planning for. Clustering of statitw pairs on any particular 
side of the earth results in highly unstable critical 
configurations, which introduce unacceptable correlations between 
the estiaated paraaeters. In this case, for instance, the 
Z-coBponent of the station positions is extreaely weakly 
deterained, being highly correlated with the T-coaponent, a 
result of the fact that all stations are at about the saae 
latitude. A coaparison of the results obtained froa 
range-differences with those obtained froa the saae data set but 
using the range approach instead, shoved that the latter is aore 
severely affected by the poor geoaetry than is the foraer. 

It is hoped that a new set of data, which is now being 
selected froa the hEfilT *80 Lageos data, will contain a 
sufficient nuaber of observations froa other areas of the world 
<e.g., the Pacific Ocean and the Australian stations) which will 
allow for a aore aeaningful coaparison of the actnal estiaates to 
those obtained by other techniques. 

Tan Gelder,B.H.w. (1978) "hstiaability and Siaple Dynaaical 
Analyses of fiange (Bange-Bate and Bange-Difference) 
Observations to Artificial Satellites," Dept, of Geodetic 
Science Bep. 284, Ohio State Oniv., Colnabus. 


30 



2.22 Doppler Experiments 


2.221 Geometric solution using ranges derived through simul- 
taneous Doppler observations: Victorlaville data set 

Introduction 

The National Geodetic Survey (N6S) was Invited by Sheltech Ltd.. 
Calgary, Alberta, Canada, and Quebec Lands and Forests to participate 
In the evaulatlon of the Doppler data set collected by eight stations 
that were observing during the period of October 2-5, 1979. These 
stations were occupied simultaneously using JMR-1 receivers. The sta- 
tions were spaced within 6 - 30 km of each other. These stations have 
been tied by conventional survey techniques to the national NAD net- 
work. The survey was done under second-order specifications. Stations 
E43.404and 273 were established on existing points. 

Doppler Test for Point Positioning. 

The above data were processed and reduced by NGS with point- 
positioning program DOPPLR and the "precise" ephemerides. These re- 
sults were used to verify relative positioning accuracies obtained 
by other Doppler reduction programs— multi-station short arc programs 
GEOOOP and SAGA III. 

The same data were also processed and reduced properly to be 
used by the GEODOR relative positioning program through a geometric 
range adjustment (geometrical solution). 

This report focuses only on the data reductions performed to be 
used by GEODOR, on the results and comparisons with terrestrial net- 
work positional information. 

Data Processing and Reduction 

The Doppler data sets were received from Richard Moreau, Direc- 
tor of Service de la Recherche et du Developpement Technologique. 

The raw data had been preprocessed by Sheltech using a program simi- 
lar to OMR's SP-2T version. 

The data that were received are: 

—one tape containing the Doppler data, Victorlaville project 
—meteorological data 
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—page from GEOOOP showing Initial values of receiver delays* approx- 
imate coordinates and antenna height 
—pages 11 and 12 from the Sheltech report 
—copy of letter to P. Gagnon giving the. final precise coordinates 
obtained by standard methods, undulations as computed by G. 
Lachapelle, and (added by hand) antenna heights. The undulations 
refer to the Clarke 1866, and the datum shifts with respect to 
1927 NAD (nominal Canadian values) are: 

AX a 15 m 
Ay « -165 m 
Az « -175 m 

The problem of positioning with geometric solution requires six 
simultaneously observing stations. In this data, eight stations at 
most were simultaneously observing. Fig. 1 is a sketch of the test 
network. 

Derivation of the Ranges 

(1) The pass header timing, the mean orbital parameters, and 
the ephemeral parameters were the main input to a program supplied to 
us by H. White. The output was the state vectors at two-minute in- 
tervals for all the passes. These state vectors refer to an earth- 
fixed geocentric coordinate system. 

Decoding the ephemeral parameters, it turned out that the data 
(as they were majority voted) have ephemeral parameter information 
for at most two minutes before lock-on time. In some passes, it was 
found that in the ephemeral parameter information the time was wrong. 
These passes were taken out. 

(2) Using the remaining passes, the Doppler counts and the 
corresponding time information was the input to the program "DOPPLRl." 
This program is a modification of a part of the proyram JMRDOP-A 
supplied to us by Larry Hothem. The output was station ID, satellite 
ei time, duration of the Doppler counts, tz time. The above informa- 
tion were sorted according to satellite and time. 
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Examining the above results* It was found that there were dis- 
crepancies In the time corresponding to the duration of the Doppler 
counts. In order to correct these discrepancies, the following pro- 
cedure was applied: 

(2a) Given the satellite coordinates at two-minute Intervals 
(taken from White's output), the polynomial coefficients were computed 
for each pass and each coordinate separately. This was done by fit- 
ting each position component as a function of t^ to a fifth-order 
polynomial. 

The general equation for X Is given by 

* ao + aiti + a2ti + ... + (I) 

The same equation was applied for the Y and Z coordinates. 

Using the fifth-order polynomial coefficient for e<«ch pass 
and each coordinate (as obtained from step (2a)), the coordinates 
of the satellite were able to be computed at the epochs ti, t 2 (see 
Fig. 2). Knowing the coordinates of the ground stations, the distances 
Si and $2 were computed. Therefore, It follows that: 

t2 ■ ti + Ti + d + ei + T - C2 - d - Tj (2) 

t2 - ti » Ti + ei + T - 62 - T2 (3) 

where 

t 2 - ti - 4.601016 s or t 2 - ti * 4.9746028 
d * receiver electronic delay 

The values ti, ei, T, ei, xz were computed using the data. So the 
computed value 

a ■ n + ei + T “ £2 - t 2 

Is different from the value t 2 - ti. Therefore, the correction b 
to be applied to the T-term Is given as: 

b ■ (ti - ti) - a 

The following assumptions are made: 

(a) No errors In the satellite clock, 

(b) No relativistic phenomena. 
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(c) The only scjrce of errors Is the duration of the Doppler counts 
T as computed from DOPPLRl (JMRDOP-'A modified program). 

These assumptions will be elaborated on In a report now In preparation. 

(3) Using the above Information, a program was written to com- 
pute the geometric ranges for each observation for all the passes. 

The Input for this program was (a) station ID, satellite ID, begin- 
ning of each observation, Doppler counts, ci time, duration of the Doppier 
counts, eatime, correction to be applied to the duration of the 

Doppler counts; (b) fifth-order polynomial coefficients for each pass 
and each coordinate of the satellite. 

The output was: station ID, satellite ID, beginning of the 
Doopler counts, Doppler counts, duration of the Doppler counts, 
correction to the duration of the Doppler counts, X,Y,Z satellite 
coordinates, and the geometric range. 

The Input to the program was sorted according to station ID 
and satellite ID. 

(4) The "OVERLAP*' program has been used and the connon obser- 
vations have been Identified. OVERLAP was fed the output of step 
(3) above, after sorting It according to satellite ID and the time. 

The output of OVERLAP was: 

(a) pass number, number of co-observing stations, satellite ID, 

(b) year-month-day, hour-minute-seconds, modified Julian day for the 
beginning of each pass, 

(c) year-month-day, hour-minute-seconds, modified Julian day for the 
end of each pass, 

(d) ID of co-observing stations, 

(e) number of events of the co-observing stations. 

(5) A program was written to select the common observations 
where six or more stations were observing simultaneously. The Input 
to this program was: 

(a) the output of (3) above after sorting It according to satellite 
and the time referred to the beginning of the Doppler counts, 
and 
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(b) the output of the fourth step (OVERLAP output). 

The output of this program was: station ID, satellite ID, modi* 
fled Julian day for the beginning of each Doppler count, Doppler 
counts, duration of the Doppler counts, correction to the duration of 
the Doppler counts, X,Y,Z coordinates of the satellite, and geometric 
range. 


(6) The meteorological data for each station were sent to us 
separately as they were recorded during the period of the observations. 

A program was written to get the meteorological data for each station 
at the mid-arc epoch of each pass (using linear interpolation). The 
input if this program was: 

(a) the OVERLAP output as it was obtained 1n the fourth step (taking 
into consideration only passes which contained six or more co- 
observed stations), 

(b) station ID, time in modified Julian day, pressure In millibars, 
temperature In degrees Celsius, relative humidity. 

The output of this program was: 

(a) the OVERLAP output as it was obtained in the fourth step (consid- 
ering only passes which contained six or more co-observed stations), 
and 

(b) station ID, modified Julian day at the mid-arc epoch of each pass, 
pressure, temperature, relative humidity referred to the mid-arc 
epoch of each pass. 

(7) A final program was written to obtain the Doppler range 
differences, the tropospheric refraction correction, and the input 
files for the 6E0D0R relative positioning program. The input of this 
program was: 

(a) the output of the sixth step, 

(b) the output of the fifth step. 

The output of this program was: 

(a) satellite ID, modified Julian day (integer), time of the beginning 
of each passj time of the end of each pass, number of stations 
observing in that particular pass, ID of the pass, number of 
events for the particular pass, maximum number of stations observ- 
ing this particular pass. 
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(b) station ID, ID of the pass, o term, b term, time from the begin- 
ning of the pass. The above Information was referred per pass 
per station. 

(c) time In modified Julian day, X,Y,Z satellite coordinates, number 
of co-observing stations per event. This Information was referred 
per each event. 

(d) time In modified Julian day. geometric derived range, elevation 
angle of line of sight to satellite, Doppler derived range, and 
ID of the station. This Information was referred per each obser- 
vation. 

The output of the seventh step was obtained by applying the following 

formulas: 

r « Ar + ro + A^^ 

where: 

Ar « \o [N^ - (fio - foo)t^] 

AO » c/foo * wavelength of adopted frequency of satellite oscillator 
c » the speed of light (299792458. m/s) 
foo “ adopted value of local reference frequency (400 MHz) 
foo * adopted value of transmitted frequency by the satellite 
(399,968 MHz) 

* time Interval between two observations 

’ Doppler counts corresponding to time 

ro » + {^ 5 “Yq)" + (Z^-Zq)^ » Initial range obtained through 

the Initial state vectors geometrically and therefore biased 
Afr * tropospheric refraction correction » (no-l)Ho/s1n E 

no » 1 + 10'® [N^ + N^] In which and denote the dry and wet 
contributions to the total refractive Index 
= 77.6 Po/To 

» 3.73 X 10® 

where Po = atmospheric pressure (millibars) 

To » temperature (degrees Kelvin) 
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eo “ water vapor pressure » 

• 6.11 X R X io7-5(T.-273.3)/(T,-35.8) 

R * relative humidity (percentage) 

Hq ■ scale height ■ 29.2 (Tq-30) 

E « elevation angle of line of sight to satellite, 
therefore changing for each observation 


(8) Using the output of the seventh step the GEOOOR relative posi- 
tioning program was used. This program was based on the geometric range 
adjustment. The mathematical model used for the range adjustment is the 
following: 

r - [(Xs-Xs^)» . * *i * “l *1 


where 

^S* ^S* earth-fixed satellite coordinates 

Xp , Yp , Zp are earth-fixed station coordinates 

Gi G.j 6^- 


a. 

1 




is a bias for each pass and station in the computation 
of the initial range. This term can take care of other 
constant biases such as satellite oscillator offset and 
receiver oscillator offset 
- time dependent bias 

= time from the beginning of each pass 


The “GEODOR" program was run 

(a) with a minimum constraint solution. The output showed that the reduced 
normal equation matrix was singular. 

(b) fixing X, Y, Z for the three stations 404, 273, E43. These three stations 
were established on existing points. The output showed again that the 
reduced normal equation matrix was singular. 

(c) fixing again the X, Y, Z for the above stations and the initial and 
final range to somehow stabilize the orbital arc. The output showed 
again that the reduced normal equation matrix was singular. 

The problem seems to be the geometrical configuration. Simulated data 
will be used to verify this guess. The procedure for the range adjustment 
is based on first-order partitioning regression, and it was an interative 
procedure described OSU Dept, of Geodetic Science Rep. 199. 
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2.222 Geometric solution using ranges derived through simultaneous 
Doppler observations: EDOC-2 data set. 

Since the progress reported on in the last Semiannual Report, work 
on this project has continued and will be reported on In a report to be 
issued in the near future. 
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2.223 


Hork k&s been coapleted on converting tke Canadian Geodetic 
Sncvey (CGS) GBOOQP pxograns for operation here at OSO. is 
described in the last seninnnoal report, this systen of prograns. 
which was written aostly bj Br. Jan Kouba of the CGS. was coded 
in CDC Fortran iv. it has been converted here to operate using 
an IBB Fortran conpiler. specifically the IBB B*Bxtended 
Snhanced. or Q conpiler (with optinization level 3 and 
lOTOBBL (OBLFiB) options in effect) • 


TkS,SBgP9£-gg9ffMB. 

The purpose of this sjstea is to process Doppler 
observations, taken with several possible types of instrunents. 
of the Mavy navigation Satellite Systea (nnss) . The nanes of the 
specific prograns in the ^stea and their overall functions 
follow: 

1. FBEDOP - Preprocesses and edits najority voted data froa 
JBS. CB1*-722. C1U-7S1. and raw data froa flX-1502 receivers. 

2. HHLFIT - Creates a (binary) file of curve*-fitted Chebyshev 
coefficients froa Haval Heapons Laboratory (MHL) precise 
satellite epheaeris. 

3. BBB6B - Berges (up to 15) single stations into one 
aulti-station file, or will serge one or nore station files with 
NBL fitted precise epheaeris. 

4. PRBPAB - Preprocesses Geoceiver or TBIHBT data, and serges 
it with either broadcast or precise satellite epheaeris provided 
by PBBDOP or MBLFIT respectively. 

5. PfiSSBO - Lists a PREOOP. BERG£. or PRBPkB created binary 
file. 

6. MiLDOB - Lists a HHLFIT created binary file. 

7. GSOOOP - Does a pass by pass sequential adjustaent of the 
preprocessed data for up to IS stations, and provides a solution 
of station positions (Z. X. Z. <). X. h). receiver inforaation. 
orbital and refraction biases, along with appropriate 
variance-covariance estiaates . 

Bach of these prograas has been converted to operate on our 
systea. Considerable testing has been perforaed to insure that 
the IBB version gives results which are identical to the CDC 
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vorsioa (with thw iaplicit assaaption that the CDC wecsion gives 
the CQccect reaolts) . Cacrently soee portioos of the systABy 
although coawerted, reeain uatested, specifically in cases where 
test data is uaavailahle aad where data of such type will not be 
used here (i.e.« nsPUl and portions of PBSOOP have not been 
tested) . The systea is cospletelj tested for JHB, 011-751 
(foraat 1H) aad BX-1S02 data. 


Description of the Conversion to IBB Fortran 

The conversion itself proceeded as follows: An initial copy 
of the systea (designated GEOOOP IV) was provided by J. Kouba in 
Decesber# I960. The majority of the conversion to IBB Fortran 
was coaplete by Ilarch, 1981* but several differences appeared 
between the IBB aad CDC version results. Br. Kouba was extreuely 
helpful here in that he made runs which allowed detailed cheching 
of the 6B00QP piograa itself. It was found that a few minor 
logical errors in the original program were basically at fault« 
as the IBB coapiler interpreted a few prograu segsents 
differently. Also, xt was finally realized that the CDC version 
sent here was not identical to the version being used by Hr. 
Kouba. Hhen a copy was sent bach to him, his runs gave identical 
results for the test data set of OSO CBA-151 data. 


In July, 1981, Br. Kouba sent his revised version of the 
6B0D0P systea, designated GBODOP V, which contained several 
corrections to the GBODOP systea, most notably in the ionospheric 
refraction correction, the partial equations, and other 
corrections that were noted here during the conversion. Several 
options had also been added including a gravity nodel, a receiver 
tiaing adjnstaent in rasDOP, and aore inportantly an option to 
edit out non-coobserved ("30 second") observations in GBODOP for 
uulti-station solutions. The iaportance of the corrections and 
the usefulness of the additional options led us to update our IBB 
version to the GBODOP V level. This updating was coupleted by 
late Angust, including testing with two standard data sets (the 
OSO CBA-7S1 data, and CGS JBB-1 data) and our BX-1502 data. 


In Septeaber, a few further a in or changes were received from 
Br. Kouba, which dealt nostly with printer foruatting changes, 
and the allowance of a higher degree gravity nodel in PBBDOP. 
These changes have not yet been included in the IBB version but 
probably will in the near future. Additionaly, a program written 
here called PLIGBO is now being tested, which will plot the 
results of the pass by pass GBODOP adiustuent on an electrostatic 
(Versatec) plotter. 
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giaal_Plfttta_for th« Proaraa System 


Once the Septeebec progrne changes have been ispleaented, 
along with any other changes found necessary during the 
processing described beloWf final docnaentation will be written 
explaining the differences in operation of the IBH version as 
cos pared to the CDC version, host of the changes in input have 
already been docanented and a set of ezasple runs have been 
created to allow testing and easier use of the progress. 
Docnaentation will also be written for the PLTGBO prograa. The 
IBB version of the 6B0D0P prograa will be sade available to Br. 
Kottba along with the above nentioned docnaentation. with Br. 
Koaba*s peraission, a nearly final copy has already been supplied 
to the Defense flapping Agency through Or. B. Kunar. 


Data Processing 

Since the conversion of the 6S0D0P Systes is now essentially 
coaplete, processing of the data obtained during our 
intercosparison test in October '>Hoveiber« 1979, has begun. The 
cosparison of the three types of instrusents, the JBB-lh, the 
CBA'‘751, and the BX-^1502 will be based on the answers to the 
following gnestions and planned G£OOOP runs: 

(1) Bow long will it take to reach a specified accuracy in 
position in terns of tine and/or satellite passes? To answer 
this, all of the data for each instmnent will be processed, 
priaarily to detenine how the final absolute position converges 
for each instruaeat. Also, when observing over the saae tine 
period, what is the final precision (variance) of the absolute 
position which is obtained (with G20DOP) 7 The results here will 
take into account that each of the receivers occupied a total of 
four stations, and assure there was no bias due to antenna 
positions. 

(2) The above questions can also be posed with the 
condition that only otaseriatioBs nade siaultaneously by all the 
instrunents (an "observation* iaplying a "30 second* Doppler 
count) should be processed. This will eUainate any bias caused 
by interference or loss of signal at one receiver, such as due to 
local interference, the antenna or pre~asplifier design, 
algorithn used to find and track a satellite, etc. This will 
allow a nore direct cosparison of the receivers* oscillators and 
electronics. This type of run can be nade using one of the (new) 
options of GBODOP itself. 

(3) Is there any dependence of the above results on the 
epheneris type? Identical runs to the above will be nade using 
the precise epheneris (for MISS satellites It and 19 only) . 
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(4) Hov stable are the tscelYex; oscillators? The resalts of 
the above rasa can be ezasined to detersine the frequency drift# 
and precision of the frequency offset. 

(5) Do the receivers pick up the sane nnnber of and quality 
of satellite passus? The nunber of passes (and/or observations) 
tracked# accepted during nnjority voting# accepted by PBSOOP# and 
accepted by 6BOOOP vill be conpared. This should indicate which 
receivers have a better tracking algorithn and/or better 
antennas. 

(6) Bow well do the receivers recover the tining 
infornation fros the satellite? It nay be possible to answer 
this question with the resalts of the (new) tining adjustaent of 
PfiSOOP. However# the way in which the JHB and B1-1S02 record 
lock'-on tines of satellite passes# and the nuserous annual 
resettings of our receiver clocks say not allow this. 

(7) How do the final absolute positions of the above runs 
coapare with the known position (s) of the stations? The conputed 
positions can be ooa pared with the astronoaic and TBAHET 
deteraiaed positions of the stations (after naking the 
appropriate geodetic corrections) • 

(B) How do the relative positions coapare? Do the 
receivers have a constant bias froa the correct absolute position 
or froa each other? Do the positions change with tine? A 
coaparison of the results froa the above runs should answer this. 
Additionally# a grand solution of all data and four stations 
could be aade# and the differences between it and the point 
positioning results noted. 

(9) Are the results dependent on the tropospheric or 
ionospheric refraction nodel# or the cutoff angle of 
observations? Several runs with different options and the sane 
data set could be aade to answer this. 


Of coarse# as each of these runs are being aade# and the 
above questions answered# aodifications will be aade in the 
planning of the runs to obtain the aost useful results. It is 
hoped that the runs (and therefore the intercoaparison) will be 
coa plated in the near future# and soae infornation obtained on 
the accuracies obtainable with each instruaent type. 
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Introduction 

The objectives of this study are the development and implementation of 
the technique of range differencing with Lageos ranges in order to obtain more 
accurate estimates of baseline lengths and polar motion variations. It is ex- 
pected, and the simulations done to date confirm this, that by means of differ- 
encing quasi -simultaneous range observations a great deal of orbital biases can 
be eliminated resulting in an estimate which is virtually bias free. In the 
past, attempts were made to use real data in order to assert the conclusions 
made on the basis of the simulations. Due to the poor geometry and distribution 
of the data, these attempts fell short of providing any conclusive results. It 
was realized, however, that even with the best possible geometry and distribu- 
tion of the observations, certain physical phenomena would have to be included 
in the model. We are currently investigating these effects on the recovered 
parameters. 

Data Preprocessing 

Since it is quite improbable, if not impossible, to obtain exactly simul- 
taneous laser observations from two ground stations to a satellite, the required 
observations must be interpolated from available commonly tracked passes. Modern 
lasers have high repetition rates and given fair weather conditions and accurate 
predictions, a sequence of ranges at a rate of about 1 pps can easily be achieved. 
Since an average Lageos pass lasts about half an hour, this implies a large 
amount of data. The high altitude of the target makes it possible to track it 



from several stations simultaneously, even If the stations' separation Is of 
continental extent. The current procedure to obtain the quasi -simultaneous 
ranges [Fig. 1] from data sets such as described above Is to determine the 
"overlap" observations for the station pairs Involved, Isolate these observations 
and determine which of the two stations In each pair has a denser data distribu- 
tion. Once this Is determined, the data of the station with the most observa- 
tions are fed Into a cubic-spline Interpolator and ranges are obtained for each 
of the data points In the alternate station's record. These ranges then are 
differenced to produce the simultaneous range-difference data. The data used In 
this procedure have already been corrected for systematic errors. Figs. 2 and 3 
depict the data record for two passes of Lageos. The bars Indicate the epochs 
when the actual observations occurred and the curve which Joins their centers Is 
the spline fit to these data. The stations marked with (*) are the ones for 
which we used the original observations In the formation of the range differences. 

Sensitivity Analysis 

Our Initial sensitivity study for the proposed method Indicated that the 
key factors In the success of the method (for polar motion variations determina- 
tion) Is the geometry of the network, length of the baselines between co-observing 
stations and their absolute location. Examining the structure of the sensitivity 
matrix for these parameters we could deduce a number of Important Items. With 
proper choice of the station pairs a nearly orthogonal system can be achieved. 

This means that each pair will be primarily sensitive to one of the components 
of the variation. This Is Illustrated In Fig. 4a & b and Fig. 5a & b where the 
matrix values have been plotted for the i and n components respectively. The 
top plots correspond to a 1000 km long north-south baseline, while the ones at 
the bottom to a baseline of the same length but In an east-west direction. The 
same conclusion can be reached looking at Fig. 6a & b and Fig. 7a & b which 
depict the situation for baselines 4000 km long. By comparing the two sets of 
plots It Is evident that there Is a linear relationship between the baseline 
length and the sensitivity of the system In the variations of the pole. The 
baselines for these graphs are all situated near the O'* longitude; similar 
graphs for these baselines have been examined for longitude 90® W. The con- 
clusions for that case are the same as the previous ones If the ^ and n com- 
ponents are reversed (see Fig. 8). 
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Fig. 1 Simultaneous laser range-differences 

(The dashed lines are the Interpolated ranges.) 
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It is quite unfortunate that although several of the existing laser 
stations fulfill the optimization requirements, commonly observed passes between 
them are very sparse. The lack of a long time period during which these stations 
co-observe Lageos makes It impossible to compare these results with those obtained 
by other techniques. 

Reference Frames 

The position of the observing stations Is defined In a geocentric earth- 
fixed system which Is materialized by a mean pole (e.g., CIO) and a mean Green- 
wich meridian (e.g., BIH). The observations are tagged with UTC epochs. The 
link between this rystem and the Inertial system Is provided by the Greenwich 
hour angle of the true vernal equinox and the coordinates of the true celestial 
pole with respect to the mean pole used. The link between the Inertial frame 
for the Integration of the orbit and the true of date frame at the observational 
epochs Is provided by the precession and nutation theories of Lleske [1979] 
and Wahr [1979] respectively. The time scales Involved, UTl and UTC are related 
to each other using the difference a(UT 1 - UTC) as obtained from JPL [Fllegel, 
1981]. 

Gravitational Perturbations 

The perturbations due to the Earth's gravitational field are computed 
based on the formulation given In [Cappellarl et al., 1976]. There Is, however, 
a subtle point In the coordinate system transformations Involved which deserves 
further explanation, especially In view of the existing confusion [Reigber, 1981]. 
This will be dealt with In the next section. The point mass accelerations due to 
the attraction of the moon, the sun, and the planets Venus, Mars, Jupiter and 
Saturn can be Included at option In our current software. The planetary ephem- 
erldes, as well as the lunar and solar position vectors, are obtained by Inter- 
polation from the JPL tape containing the DE114 Solar. Lunar and Planetary 
Ephemerls [Standi sh, 1981]. 

Computation of the Non-Spherical Effects from the Terrestrial Gravitation 

The gravitational potential Is expressed as a series expansion in 
spherical harmonics. With the 0^^ -order term (point mass effect) removed, 
we can write the, non-spherical part of it as: 
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v(r.*.x) . fl V [S^sln "IX + C^cos mA] P, 


mn 


(sln^) 


( 1 ) 


where r»^,X are the spherical coordinates of the point of evaluation. From 
the purely mathematical point of view the choice or definition of the (r,a*v) 
system 1s Irrelevant. The function V, the potential » which Is approximated 
by the series is Invariant with respect to any similarity transformation of 
the underlying coordinate system; the value of V at P(r,^,A), Vp will be the 
same whether P Is defined In the (r.^.x) system or. say, the (r‘,^',.\') 
system. The series, however. In (r'.^'.x') will have different constants, 
harmonics. In this case. It is thus obvious that the similarity transforma- 
tion between (r,^,x) and (r',^*,x') propagates as a similarity transformation 


between (C^, and (C' 
nm nm •* 


S' ) 
nm* ^nm^‘ 


So far, the above Is just a restatement of already known facts. For 
Instance, Kleusberg [1980] has developed the (rather cumbersome) formulae 
to obtain the primed harmonics from the original, given the transformation 
parameters. It Is the last statement In the previous paragraph that causes 
the problems: The fact that a change In the coefficients cannot be attributed 
to an actual change of the coefficient or a coordinate system change, unless 
It Is known a priori which of the two happeneo. In a sense this Is similar 
to "the principle of equivalence" In general relativity, whereby gravitational 
and Inertial forces are Inseparable. Considering now that In a dynamical solu- 
tion the coordinate system definition Is provided by the satellite dynamics 
(short of a longitude definition), one realizes how Important It Is to clearly 
define a priori the system In which the harmonics of the series In (1) are 
referenced. 

It Is a well known face [Helskanen and Moritz, 1967; Nagel, 1976] that 
for the low-degree harmonics (up to n«2) we can easily associate with them 
certain properties of the estimated model. The pair (C 2 j, Sjj) Is of partic- 
ular Interest to us since their values are directly proportional to the align- 
ment of the reference system with the principal axis of maximum moment of 
Inertia. Lambeck first called attention to the Implications of this otherwise 
Innocent looking property. In [Lambeck, 1971], where he stated that 


u* ^ ^ 


and 


t 


( 2 ) 


u* and V* being the coordinates of the axis of figure with respect to the 
third axis of the reference coordinate system (the axes u*, v* are defined 
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In the same sense as the for the pole). A proof for (2) can be found 

in [Nagel, 1976]. This interesting property becomes a real problem in practice, 
due to the fact that the axis of figure is time variant. In [McClure, 1973; 
Nagel, 1976; Leick, 1978; and Moritz, 1979], one will find detailed descriptions 
of the motions of this axis. It suffices here to mention that there is a free 
motion with an amplitude of 2 m and a period equal to that of the Chandlerian 
wobble (“ASO days) and a forced diurnal motion with an amplitude equal to 60 mi 
To include at least the long period effect in the harmonics and 
$ 2 ^ implies that the coordinates of the center of the Chandlerian wobble are 
known. Since there are only known from observations, one has to rely on some 
predictions based on the past motion of the center. Given those, then the 
instantaneous values of the two harmonics can be computed from relatively 
single formulae: 


cJl “ lx (Xp - x)J Cjo 

» - i? ♦ I (yp - 7» Cjo 


(3) 


From some simple computations, the effect of the aforementioned motion 
of the axis of figure on 62 ^ and S 21 seems to be capable enough to cause apparent 
accelerations of the order of 10 “^^ m/sec^ on Lageos, for a change of a few tens 
of milli-arcseconus in the polar coordinates (a value quite reasonable for a 
five-day interval). We are currently investigating the implications of such an 
acceleration. 


Results and Conclusion 

As it was mentioned before, although very realistic simulation studies 
have indicated that the proposed method can produce virtually bias free polar 
motion variations, the lack of well distributed real data makes it impossible 
to deduce anything conclusive about it. As far as formal precisions are con- 
cerned, it seems quite definite by now that a few milli-arcseconds can easily 
be achieved even for one-day averages. That, of course, is only a lower bound 
of no great significance as far as accuracy is concerned, but if biases can be 
minimized with this method then the rtcl accuracy level will not be exceedingly 


U 


higher than this estimate. We have lately acquired the laser data collected 
during the short MERIT campaign, and it seems as if there are a number of more 
suitable data sets which can be used for testing the method, than that which 
we had earlier from the late 1979 to early 1980 data. Once this study is com- 
pleted, we should have more definitive answers to present. 
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3 


TIME 

MORNING 

8.30 

8.40 

9.00 


10.00 


TUESDAY 1 SEPTEMBER 
MORNING, 8.30 - 12.30 
SESSION II 


OPTIMAL DESIGN OF GEODETIC NETWORKS 
(ACCURACY, RELIABILITY, COSTS, ETC.) 


CHAIRMAN ; A.«C. tikAt 

(U^UU) 

REPORT OF SSG 1-59: COMPUTER AIDED DESIGN OF 
GEODETIC NETWORKS 

REVIEW PAPER : OPTIMAL DESIGN OF GEODETIC NETWORKS 


11. 1. OPTIMIZATION OF OBSERVING LOGISTICS IN GEODETIC NETWORKS 

(C^tk4>d^) 

11. 2. OPTIMAL DESIGN OF GEODETIC NETWORKS 

B^^4md4>^dt - diH^^4 

Ad^4 ( A^ 4. tU .i>m^) 

11. 3. ON THE DESIGN OF GEODETIC NETWORKS USING ITERATIVE METHODS 

X.«4. tK4»44l A. 8. M^f^4Cd*U.C 

(£Mf.d4k»d) 

11. 4. SECOND ORDER DESIGN OF GEODETIC NETWORKS 
PROBLEMS AND EXAMPLES 

0, k%dd44>dt k mdABAfAdi-dA (B^d» il*|k. ^ 


Coffee break 



4 






10.30 I I. 5. AN INVESTIGATION ON THE OPTIMIZATION OF THE SPACE 

OBJECTS OBSERVATIONS FOR THE EAST-EU.' ^EAN SATELLITE- 
TR I ANGULATION 

C Bvkt^cu>yi,ak) 


11. 6. OPTIMIZATION OF GEODETIC NETWORK 

L.ti. (Titfi, 

11. 7. OPTIMIZATION OF THE CONFIGURATION OF GEODETIC NETWORKS 
iL^Z. d^'Ck, UM.i,it4,K'A^ta4> 


II. S. INTERACTIVE NETWORK ANALYSIS 

It. A. tU.f^Lami £.}. 0.4 

(CoMAdAj 

II. 9. SOME CONSIDERATIONS ON THE OPTIMAL DESIGN OP GEODETIC 
NETWORKS 

do.lM» (tAd» 


II. 10. SECOND ORDER DESIGN OP GEODETIC NETWORKS BY AN 

ITERATIVE APPROXIMATION OF A GIVEN CRITERION MATRIX 
a. %AcJitU.A<XA UAA-AAA'iAtit HHacLaa a^d. 

Za^, 0»4 ^A\>mAAI^) 

OTHER CONTRIBUTIONS (ONLY FOR SYMPOSIUM PROCEEDINGS) 

1. OPTIMIZATION OF A LOCAL CONTROL NET FOR GEODYNAMICS 

BAAdAcJuAAl U., UAt*^ck, OAt^t^ckA^ <fAAddtAA<M.A^ Aku^^AAit-utt^t 

HUHAAkAA ChAd, ZaA‘ A.4 ^AmAAIf) 

2. POSSIBILITIES FOR THE STRENGTH IMPROVEMENT OF THE HELLENIC FIRST 
ORDER TRIANGULATION NETWORK 

A, 0A\mtkAAA /i4. "hAtAAMkl B . iA44*M,Af*AAdAA ,U»AAAAAAAt\f A^ ZkA 4 * AtOrAAkif (^AAAAA 

3. EXPERIENCES IN APPLYING GEODETIC NETWORK OPTIMIZATION 

Z. KaaLav , ln.4i,. %A (^AdA%A-^A - BaA^AaA, (J^AitAAAA) 

4. REPORT OF SSG 4.71; OPTIMIZATION OF GEODETIC NETWORKS (PRESENTATION 
IN SSG 4.71 - MEETING) 

^.ickmAtt, UAAAAA^Atdt ZAAt^AAkA (ItAd. Za^. a{. I^AAmAAl^) 

5. TREE-SEARCH ALGORITHM FOR THE TRIANGULATION-TRILATERATION COST 
MINIMIZATION 

C, ZiO^AAA , lUkAAAAAAkt^ A^ TJkA^’iAtAtAik.A (^AAAa) 


12.30 Lunch break 
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TUESDAY 1 SEPTEMBER 
AFTERNOON , 14.30 - 17.30 

SESSION III 


I aMODERN OBSERVATION TECHNIQUES FOR TERRESTRIAL NETWORKS 

i (SATELLITE AND INERTIAL TECHNIQUES, VLBI, ALTIMETRY, ETC.) 
1 


TIME 

AFTERNOON 

14.30 CHAIRMAN ; 2.2. (Jkla, U(U>€. - 

(US>A) 

14.40 REVIEW PAPER ; SPACE TECHNIQUES FOR TERRESTRIAL NETWORKS 

(UU) 

IS. 00 III. 1. THE ATMOSPHERIC EFFECTS ON ELECTROMAGiMETIC DISTANCE 
MEASUREMENTS IN GEODETIC NETWORKS 
S. IL. o>4 

111. 2. GEODESY AND THE GLOBAL POSITIONING SYSTEM 

2.2. (US-A) 

111. 3. USE OF SATELLITE-DOPPLER OBSERVATIONS FOR GEODETIC 
NETWORKS 

U, i. ^^ 44 , a>/> Ko^t><>»4XQM 

111. 4. NEW METHODS FOR THE REDUCTION OF SATELLITE DATA 
APPLICABLE TO GEODESY 

0a|»4. a/ ^«4aa^ca4 ^a^^aaavaa^. - 

Saa P.aa^a 


16.00 Coffee break 
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16.30 III. 5. CONSEQUENCES OF GRAVSAT AND GPS: NEW CONCEPTS OF 
GEODETIC NETWORKS 

d. S-t^ck , 

Oa,^>m^ta>dLt Ch^dL. o>/ 

III. 6. ON THE USE OF ORBITAL METHODS FOR DEVELOPMENT OF 
SATELLITE GEODETIC NETWORKS 


OTHER CONTRIBUTIONS ( ONLY FOR SYMPOSIUM PROCEEDINGS) ; 

1. ESTABLISHMENT OF A BASE LINE OVER THE SEA 

r^v«.«.ce>^ 

2. REPORT ON THE WORK OF SSG 1.42: ELECTROMAGNETIC WAVE PROPAGATION 
AND REFRACTION IN THE ATMOSPHERE 

1. Jl. Ktm S-^u>tk 

(4t^4t^At'UA) ( -itk tk^ ) 

3. REPORT OF THE SSG 1.26: CONTRIBUTIONS FROM SATELLITE GEODESY 
TO TERRESTRIAL GEOMETRIC GEODESY 

2 . JiaM,ku,\.C, (^ct^tMdL) 

4. ON THE TRIGONOMETRIC LEVELLING 

Z. , l-Af¥kit^<L^k (^i-n^keutkdL) 


17,30 END OF WORKING SESSION, CONTINUATION FRIDAY SEPTEMBER 4 
IN SESSION VIII 


19.00 Visit of the Bavarian State Collection of Painting 
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WEDNESDAY 2 SEPTEMBER 

MORNING AND AFTERNOON, 8.30 - 12.30 and 14.30 - 17.00 

SESSION IV 


NETWORK ANALYSIS MODELS (BLUNDERS, 

SYSTEMATIC ERRORS, STATISTICAL 
TESTS, DEFORMATIONS, ETC.) , 

TIME 

MORNING 

8.30 CHAIRMAN t 

8.40 REVIEW PAPER ; DIFFERENT ASPECTS FOR THE ANALYSIS OF 

GEODETIC NETWORKS 

9.00 IV, 1. A REVIEW OF MODEL CHECKS AND RELIABILITY 

IV. 2. EXPERIENCES WITH A NONSTATISTICAL METHOD OF DETECTING 

OUTLIERS 

IV. 3. ADJUSTMENT BY MINIMIZING THE SUM OF ABSOLUTE RESIDUALS 

H. ^ck9, 

IVC4. RELIABILITY AND GROSS ERROR DETECTION IN PHOTCGRAMMETRIC 
BLOCKS 

%, AckA.^>mck»*^ , Ui9*,u€>^,i'otdt (^vd, 9,^ 


10.00 Coffee break 

a 



10.30 IV. 5. ON THE RELIABILITY OF 3RD AND 4TH ORDER NETWORK 

DENSIFICATION 

IV. 6. ACCURACY ANALYSIS OF THE FINNISH LASER GEODIMETER 
TRAVERSE 

y. )LaM,Lvk*»,i0 1 X. P.«kKM, 'hi»»^¥ki.'AL • 

U, S>,4.4 o>^ 

( ) 

IV. 7. COMBINED LEAST SQUARES SOLUTION OF TERRESTRIAL AND 
DOPPLER OBSERVATIONS 

IV. 8. ERROR PROPAGATION IN LEVELLING NETWORKS 

IV. 9. THE BINOMIAL DISTRIBUTION AS PROBABILITY MODEL OF 

THE DISTRIBUTION OF ERRORS OF GEODETIC OBSERVATIONS 
Ik, a. ^lAckvK-tlka, , 

(P.<^ta,*^i) 

IV. 10. A METHOD FOR DETECTING VERTICAL SOIL MOVEMENTS FROM 
SCATTERED LEVELLING NETS. APPLICATION ON FIRST 
PRECISE LEVELLING NETS OF HAMBURG 
k. X. £M.o^ka\J. , 

Ch€>i. kt>^. 

IV. 11. ON THE OPTIMIZATION OF LEVELLING NETWORKS WITH RESPECT 
TO THE DETERMINATION OF CRUSTAL MOVEMENTS 
It. , Unk4»iAV^^i»t&t MA»u^it€,K f ( 

IV. 12. PRECISE LEVELLING ACROSS ACTIVE FAULTS IN CALIFORNIA 
A.<^, , Ot.^, Xc^ei«c<.4 , 

“ S,04U,a, (HS.A) 

12.30 Lunch break 


o 


AFTERNOON 


14.30 IV. 13. THE EFFECT OF THE VARIATION OF THE DESIGN IN A LARGE 

NETWORK, DEMONSTRATED FOR BLOCK D OF RETRIG 

IV. 14. DESCRIPTION OF HORIZONTAL STRAINS AND SOME GEODETIC 
ASPECTS OF THEIR ANALYSIS 
4 (. H€>i4cL, MAcL4cku4>*. d*kt^€>4mt>kK> - 

IV. 15. LEAST SQUARES PREDICTION OF HORIZONTAL CO^.'RDINATE 
DISTORTIONS IN CANADA 

A» , S>kt,tt^cL - 

( 

IV. 16. CRUSTAL DEFORMATION OF THE AUSTRALIAN PLATE MEASURED 
BY SATELLITE LASER RANGING: PRELIMINARY RESULTS 
A. Au'^i/k^ • I £,^. Jl(a.dt>«.4>4 / 3. At>^^4ck, 

o>f( S^u-tk ( Aw‘^t,Ko,t<,Ci,) 

15.30 Coffee break 

IV. 17. THREE-DIMENSIONAL KINEMATICS OF EARTH DEFORMATION 
FROM GEODETIC OBSERVATIONS 

H.l, jl«.a.(tX^4‘C64 0*.^. ^4- 

omA- ’hf^d(/k^t\‘CeU- )U.^ca-\><.k - (Hem tecLta^d) 

IV. 13. ANALYSE STATISTIQUE DE LA DEFORMATION DES RESEAUX 
A. It. KAteemckt - 

lAemJJk ( Hxemee ) 

IV. 19. GEODETIC PREDICTION OF ACTUAL CRUSTAL DEFORMATIONS 
AT THE SEISMIC AREA OF VOLVI 
A. Oe^mmmde/ C. I 0, ^o>4 4 / 

O', iLkmcki^e ^ XA<4>44<i»i<n*«<i*A>c (^^eeee) 

IV. 20. STRENGTH ANALYSIS OF ANGULAR ANBLOCK NETWORKS WITH 
DISTANCE MEASUREMENTS ALONG THE PERIMETER 

P,. ked^et, C^Km% (Amet^i,^) 
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OTHER CONTRIBUTIONS (ONLY FOR SYMPOSIUM PROCEEDINGS) 


1. ON THE POINT IDENTITY PROBLEM IN CONTROL NETWORKS 

a.. Cha,d. 

2. TECHNIQUES OF STATISTICAL ANALYSIS AS APPLIED TO GEODETIC 
ADJUSTMENT COMPUTATIONS 

H,k. JIa4€AAaL & OAAAtAfUHAAk Ot>V , - 

Hj^d,A\>AA>Ad> (lAd-OA) 

3. HEIGHT DEPENDENT ERRORS IN SOUTHERN CALIFORNIA LEVELING DATA 

0,0. ^AcL 4A», Ot.^, tAAtk AaJL ici.AACA6 “ £a* (UiA) 

4. REPORT OF SSG 4,60: STATISTICAL METHODS FOR ESTIMATION AND 
TESTING OF GEODETIC DATA 

( P,AA6AIAtAl.AAA AA tkA S.S-^“PtAAtAA^ ) 

H.Z, ILaaL, tlAAAAAAAt&k 3 AAA (^AiL, RaIi^ . a4- ^A\MAA>^) 

5. OPTIMIERUNG DES MESSPROZESSES 
(OPTIMIZATION OF THE MEASUREMENT PROCESSES) 

R, iRijAAAt %A ^AaJLaaA^A - BaA^AaJ, (t^U^AAtAAAA) 

6. STRAIN OF HORIZONTAL NETWORKS 

R, UaaACAL, Oa^U.. A^ iAAAA^AA^ • %AAiAAACtAA (CaaACLa) 


17.30 City Tour 


TIME 

MORNING 

8.30 

8.40 

9.00 


10.00 


THURSDAY 3 SEPTEMBER 
MORNING, 8.30 - 12.30 
SESSION V 


ADJUSTMENT PROCEDURES OF GEODETIC NETWORKS 
(POINT AND INTERVAL ESTIMATION, VARIANCE AND 
COVARIANCE COMPONENT ESTIMATION, ETC.) 


CHAIRMAN ; C. 2. • 

UckiiUU (UU) 

REVIEW PAPER: ADJUSTMENT PROCEDURES OF GEODETIC NETWORKS 
('Sc.i. a/. ^KmAt^f.) 

V.l. MATHEMATICAL MODEL FOR NETWORK GEOMETRIC ADJUSTMENT 

C. tt. Oac^/ C. Aa'6«>a»a4 - 

V.2. THE ESTIMATION OF COVARIANCE MATRICES FOR PHOTOGRAMMETRIC 
IMAGE COORDINATES 
a. Ji. ^cLK4>t,L, 

£a(». a4 ^AAAAA^^ 

V.3. ADJUSTMENT OF GEODETIC NETWORKS IN SPACE 

4 , ikAAA44A, a/ Xa^4^ 

V.4. ON FREE NET ADJUSTMENT ON THE ELLIPSOID 

iK. XcXh^A^IacX, IaC^AA4a4a i(^A«AA4A4A4 
Ch*^. Z€>t^, A^ ^AAMAAykJ 


Coffee break 
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10.30 V.5. ON THE ADJUSTMENT OF CONTINENTAL VERTICAL GEOKINETIC 

NETWORKS 

2.)t. - 

) 

V.6.. COMBINED ADJUSTMENT OF DOPPLER AND TERRESTRIAL NETS 
IN DOPPLER REFERENCE SYSTEM 

2 . ^40^4,4^ ?.a4^^^aMm4,t4,f 

AJIZ - 0t4%t^t^ 

V.7. AN APPLICATION OF THE GRAPH THEORY IN THE ADJUSTMENT 
OF A GEODETIC NET 

Z, • Ci,4>44,tt4,i , I4t4^i>t4t4> 4> 

- UtkiuA^,4^tA Z^4,4,4t4, 

V.a. JUNCTION OF EUROPEAN DOPPLER OBSERVATION CAMPAIGNS 
CARRIED OUT BY THE IfAG/SFB 78 

V.9. AN INTRODUCTION TO THE ADJUSTMENT OF ASTRO-GEODETIC 
NETWORK OF THE PEOPLE'S REPUBLIC OF CHINA 

0. ^(4, C/i^t44’44>4C^A^ (^44,tt4CkA^ ^ 2 . 4 . UkUtJL 

\L4>tli44^4^^»ktiJlU CJH 0 I 44 , - BA'LmM.%kl4At4^, 

( Cki0t44>) 

OTHER CONTRIBUTIONS (ONLY FOR SYMPOSIUM PRQCFSDTNGS ^ : 

1. REPORT OF SSG 1.69: COMPARISON OF CONTROL NETWORK ADJUSTMENT MODELS 

( 4 . A4kU.*.li4A%<0 , Ut4^iU.\4^f‘ kc^tt'CM^kAm ( 

2. READJUSTMENT OF THE POLISH PRIMARY HORIZONTAL NETWORK 

^%i%40<.iu<0 , C^t^kKtum ’hl4^Km44>kf:€0%t44t ^ - 

HA^4Am (P^tojuL) 

3. ANWENDUNG DER KOLLOKATION fOR ZWEI GEMEINSAM ABHANGIGE STUTZ^/ERTE 
UNO FUR DIE IN DER EBENE LIEGENDEN STOtZPUNKTE (ANWENDUNG DER 
KOLLOKATION ZUR BESTIMMUNG DER OBJEKTIVDISTORTION) 

)•. , Z4,L\>4t>t4kt fL. M&Lt,K4> <^^A&440t0^P,K4M,Ck i t%4.<,k4f4t^^i.aJt0i,4.) 


12.30 Lunch brea»c 
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THURSDAY 3 SEPTEMBER 
AFTERNOON, 14.30 - 17,30 
SESSION VI 


COMBINATION OF HORIZONTAL, VERTICAL AND GRAVITY NEP/ifCRKS 
(THREE AND FOUR DIMENSIONAL PROBLEMS, REDUCTIONS , ETC. ) 


TIME 

AFTERNOON 

14.30 CHAIRMAN ; 

( ) 

14.40 REVIEW PAPER ; COMBINATION OF HORIZONTAL, VERTICAL 

AND GRAVITY NETWORKS - A REVIEW 


15.00 VI. 1. COMBINATION OF LEVELLING AND GRAVITY DATA FOR 

DETECTING REAL CRUSTAL MOVEMENTS 
B, lU^t4K^L4, 

VI. 2. A CONTRIBUTION TO 3D-OPERATIONAL GEODESY 

f>.d, iU>^, Ch^dL, 


VI. 3. ORIENTATION INFORMATION OF LEVELLING AND GRAVITY 

MEASUREMENTS IN THREE-DIMENSIONAL REGIONAL NETWORKS 
A. A^tm, D^Ut4<Ju,4 ^4,A&t,4,4<,L4.4 %4^4<,ktM^4'OI^4td»t<U. 

- (‘ht.JL* At-i^. 4>4> 


VI. 4. AN ASTRO-GRAVIMETRIC COMPUTATION OF THE QUASI- 

GEOID OF THE FEDERAL REPUBLIC OF GERMANY 

0. I 0. CkAt,4i4>l A, AomcL , 

At^^mo0t^dt>4> ^4>dA4<CA • At,^, 


16.00 Coffee break 
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16.30 


VI. 5. ON THE INTERPOLATION OF GRAVITY ANOMALIES AND 
DEFLECTIONS OF THE VERTICAL IN MOUNTAINOUS 
TERRAIN 

VI. 6. test-computatiOns of three-dimensional geodetic 

NETWORKS WITH OBSERVABLES IN GEOMETRY AND GRAVITY 

space 

Chd. 

VI. 7. THREE-DIMENSIONAL ADJUSTMENT OF GEODETIC NETWORKS 
USING GRAVITY FIELD DATA 
A, 2. kAAtt^f, ^e.a^A^4-c-c4 

S^AAAAtA^AA AaA 2lAdAA^AAAt jlAiAAAAk - AAttAA^t.AA 
(flAm ZAAtAAd) 

17,30 End of Working Session 


17.30 


20.15 


Departure to Castle Lustheim 
Serenade at Castle Schleiflheim 


FRIDAY 4 SEPTEMBER 
MORNING, 8.30 - 12.30 
SESSION VII 


COMPUTATIONAL PROBLEMS IN CLASSICAL AND NON-CLASSICAL 
ADJUSTMENT MODELS (LARGE NETWORKS, ILL CONDITIONING, 
COLLOCATION, SOFTWARE PROBLEMS, ETC.) 


TIME 

MORNING 

8.30 CHAIRMAN ; A. 

REPORT OF SSG 1.2 1:NUMERICAL COMPUTATION OF LARGE 
TRIANGULATION NETWORKS (together with A. PXxttM,) 

8.40 REVIEW PAPER ; COMPUTATIONAL PROBLEMS IN CLASSICAL AND 

NON-CLASSICAL ADJUSTMENT MODELS 

( OvuumokKU,) 


9.00 VII. 1. REPORT ON CONCLUSIONS OF THE SYMPOSIUM "MANAGEMENT OF 
GEODETIC DATA" 

r 

VII. 2. AUTOMATED SET-UP AND ADJUSTMENT OF TRAVERSE NETWORKS 
h. 3 a^^*.JLm€., X«.c4.»V4c4«. ( Au,it\>-LA) 

VII. 3. UNIVERSAL PROGRAM FOR ADJUSTMENT OF ANY GEODETIC NETWORK AND 
FOR CRUSTAL STRAINS 

X. Hcu^a-dLi^ , ^o>C'4/4c At-\^ Co>. , Ztd>. 


VII. 4. PROG.iAM SYSTEM "NET2" s'OR THE ADJUSTMENT AND ANALYSIS OF 
LARGE GEODETIC NETWORKS USING SPARSE ALGORITHMS 

(^A. a 4> ^AKmAAIf.) 


10.00 Coffee break 
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10.30 VII. 5. iFARiE MATRIX ALGORITHMS "NESTED DISSECTION" AND • 

"MINIMUM DEGREE ORDERING" APPLIED TO DHM GENERATION 


VII. 6. SOME NEW PROCEDURES OF THE SEQUENTIAL ADJUSTMENT 

A. i. ^40i^VA4iu.'6V‘4-u - 

(?.^ta,u.eL) 

VII. 7. CHOIX DE METHODE DE CALCUL EN ESTIMATION DES PARAMETRES 
SELON BAYES 

2. X'VAM.C^u/ h. yto.'V - 


VII. a. THE INCOMPLETE CHOLESKY CONJUGATE GRADIENT METHOD 
FOR NETWORK ADJUSTMENT 

Ti, 0. <Li, 

VII. 9. THE OPTIMIZATION OF LARGE TRAVERSE NETWORKS 

A.. ?j€aU.€.lc , 0>/- 0M^<L ~ 

C\>a-co»m (Po^tM-dL) 


VII. 10. COMPUTATIONAL PROBLEM IN THE INDONESIAN-MALAYSIAN 
INTERNATIONAL BOUNDARY PROJECT 

0, X. S-Af^ta-d.€.m<^, Kikt, a«^ 


VII. 11. ORDERING AND DISSECTION OF GEODETIC LEAST SQUARES 
EQUATIONS 

}L t A» A* f ”* 

( 0cKma.>Ut,) 


OTHER CONTRIBUTIONS (ONLY FOR SYMPOSIUM PROCEEDINGS) 

1. AN AUTOMATED THREE DIMENSIONAL ANALYSIS OF DEFORMATIONS FOR 
CONTROL OF BUILDINGS 

J,. 3a,Lu^<ia>^^l ' Z. (^ceL. Pc^. o<4 
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2. AN ALGORITHM FOR THE INVERSION OF THE CORRELATED OBSERVATIONS 
COVARIANCE MATRICES 

ha,c<.Mcui> dt, C<,ui,t • j!c4^a>o> 

3. A CONTRIBUTION TO THE TREATMENT OF DEFECTS IN LARGE GEODETIC 
NETWORKS 

4. LEAST SQUARES ADJUSTMENT OF LARGE-SCALE GEODETIC NETWORKS BY 
SPARSE ORTHOGONAL DECOMPOSITION 

a>^ (ILS^A) 

5. ESTIMATION DE LA CONDITIONNALITE NUMERIQUE DES MATRICES NORMALES, 
OBTENUES EN COMPENSATION DES RESEAUX DE NIVELLEMENT 


12.30 Lunch break 


FRIDAY 4 SEPTEMBER 
AFTERNOON, 14.30 - 17.00 
SESSION VIII 


CONTINUATION OF TOPIC 5 

' I 

MODERN OBSERVATION TECHNIQUES FOR TERRESTRIAL ! 
METa/ORKS (INERTIAL TECHNIQUES, VLSI, ETC.) I 


TIME 

AFTERNOON 


14.30 Vlil.l. USE OF INERTIAL TECHNIQUES FOR GEODETIC NETWORKS 

4. - (^i,mKCk(L«k (USA) 


VIII. 2. ADJUSTMENT PROBLEMS IN INERTIAL POSITIONING 

VIII. 3. WHAT INFORMATION CAN YOU GET OUT OF AN INERTIAL 
SURVEY SYSTEM, AND WHAT CAN YOU DO WITH IT? 

VIII. 4. ORIENTATION OF SETS OF DIRECTIONS BY GYROSCOPIC 
MEASUREMENTS 

Id. U, Ut.i.4t4.^/ P.. 


15.30 Coffee break 

16.00 VIII. 5. USE OF INTERFEROMETRIC METHODS FOR GEODETIC HORIZONTAL 

NETWORKS 

Pt.Z. • Ma>ACAm 

(USSZ) 

VIII. 6. MEASURES FOR PRECISION AND RELIABILITY IN PLANNING, 
ADJUSTING AND TESTING A EUROPEAN VLoI-CAMPAIGN ; 

SURVEY OF ONGOING AND FUTURE MEASUREMENTS 
k. - 0€.Li.t 
( ) 

17.00 CLOSING SESSION 
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lAG Sifinpo4-im on G&odtUc Nttuo-iki and Compototutu 
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GEODESY AND THE GLOBAL POSITIONING SYSTEM 


Ivan I. Mueller and Brent Archinal 
Oeoartment of Geodetic Science and Surveyinq 
The Ohio State Univesity 
ColumbuSt Ohio 43210 USA 


ABSTRACT. The 18 satellites at altitudes of 20 thousand kilometers proposed for the NAVSTAR 
Global Positioning System will provide an opportunity to determine the relative positions of sites 
hundreds of kilometers apart to accuracies on the order of centimeters after a few hours of obser- 
vations, This presentation reviews recent efforts In the United States In the development of dif- 
ferent geodetic receivers designed with the above goal In mind. 

1. INTRODUCTION 

This presentation, prepared at the request of the organizers of this symposium, is a review 
of the current (early sutnner, 1981) status of the NAVSTAR Global Positioning System (GPS). The 
18 satellites at heights of some 20,000 kilometers proposed for the GPS will provide an opportu- 
nity to determine relative positions of sites up to hundreds of kilometers apart to centimeter- 
type accuracy after a few hours of observation, possibly less. The GPS and its navigational use 
Is described in [Institute of Navigation, 19c0], Geodetic receivers to exploit the GPS signals 
are being developed mainly under the cognizance of a U.S. interagency coordinating group composed 
of representatives of the National Aeronautics and Space Administration (NASA), the Defense Map- 
ping Agency (DMA), the Geological Survey, and the '.ational Geodetic Survey (NGS). The most prom- 
ising geodetic receiver developments and tests are being conducted with the support of this inter- 
agency group by the Jet Propulsion Laboratory (under NASA sponsorship) [MacOoran, 1981], and by 
the C.S. Draper Laboratory (CSDL) and the Massachusetts Institute of Technology (MIT) [Counsel- 
man, 1981], sponsored jointly by the other three agencies. DM.A also continues tfO support activi- 
ties at the Naval Surface Weapons Center (NSWC) with equipment by Stanford Teleconmunications, 

Inc. [Anderlo, 1980]. Receivers are also being developed by other commercial companies on their 
own. Most no':ab1e are those by the Canadian Marconi Company, Magnavox Advanced Products and Sys- 
tems Co., and Texas Instruments, Inc. (see Table 3). Most of these developments are described 
on a continuing basis in the [CSTG Bulletins]. 

2. DESCRIPTION OF. THE GPS 

The GPS has been designed to provide for passive, all weather navigation from virtually any 
point on earth at any time. The GPS satellites will be operated by the U.S. Air Force for the 
U.S. Department of Defense. The GPS constellation of satellites (known as the "Space System Seg- 
ment") was originally to have consisted of 24 satellites equally spaced in three orbit planes, 
with 12-hour periods. This would have provided for cc'.tinuous all weather precise navigation, 
with four satellites always visible in good geometry (above five degrees altitude) for any point 
on earth. However, for budgetary reasons, the constellation has been cut to 16 satellites. 

These satellites will probably still be placed in the same orbits as originally planned, except 
that the six satellites in each orbit plane will not be equally spaced. Such a configuration 
will provide the best orbital geometry and visibility over the longest periods of time worldwide 
(although with losses of up to 40 minutes over an eight-hour period) [Jorgensen, 1980]. 


z 


The GPS satellites themselves will employ precise atomic frequency standards to assure high 
oscillator stability and precise timing information. Each satellite will broadcast on two 
L-band frequencies, (1575.4 MHz) and L2 (1227.6 MHz), with the two frequencies allowing the 
elimination of first-order Ionospheric refraction by receivers. Three different codes will be 
modulated onto the L-band frequencies, these Including a D code consisting of a satellite message 
(at the rate of 50 bits/second) and two pseudo-random noise (PRN) codes called the coarse/acqul- 
sltlon (C/A) and precision (P) codes. These latter two codes are respectively modulated onto 
the L-bands at frequencies of 1.022 MHz and 10.23 MHz, with the codes repeating every one milli- 
second for (C/A) and every 38 weeks for (P) (In practice, P will be reset every week). These 
codes serve three purposes: (a) to Identify each satellite uniquely, (b) to provide for measure- 
ment of the signal travel time to the receiver by measuring the phase shift required to correlate 
the codes to receiver generated codes [Mil liken and Zoller, I960], and (c) to provide for restric- 
tion of the use of the system If desired, since the P code must be known a priori by the user to 
obtain the highest accuracy (except for certain types of receivers described later) [Stansell, 
1980]. 

A normal GPS receiver would thereby have to complete the following operations to successful- 
ly track a satellite: {&) acquire and recognize the C/A code, (b) obtain the necessary satellite 
message Information to obtain the P code epoch by decoding the 0 code from the C/A code, (c) 
acquire and re;,wgn1ze the P code. 

The present status of the GPS is that six NAVSTAR satellites are presently In orbit (al- 
though one of them has a faulty clock). A seventh satellite Is to be launched this sunmer, and 
the entire system is to be In operation by 1987. These satellites lie In two planes Inclined at 
63 degrees to the equator, allowing four satellite observations to be made (over three hours) 
from several locations, such as the U.S., a large portion of South America, Australia, the Near 
East, and above latitudes of t 60‘ [Perreault, I960]. Therefore, although limited somewhat In 
location and time, the system Is now operational, so that testing or limited actual use of any 
GPS receivers Is a possibility. 

This description would not be complete without some mention of the ground facilities, which 
•are known as the "Control Segment" of GPS. Four ground stations will continuously monitor all 
GPS satellites In view. The master control station will be located at Vandenberg Air Force Base, 
where all tracking data will be collected, with precise orbital Information (predictions) then 
being made. This Information will then be loaded Into the satellite memory at least once per day 
as the broadcast message for the satellite. Further Information on the Control Segment Is avail- 
able In [Russell and Schalbly, 1980] which also discusses the expected accuracy of the orbital 
predictions. 

3. GPS OBSERVABLES 

Three major methods of using the GPS satellites for precise positioning .have been proposed, 
making use of pseudo-range, Doppler, and Interferometric observations. Instruments either have 
been or are being const: \>cted to observe the satellites using variations of these methods. De- 
tails of these proposals can be found In [Fell, 1980] along with extensive simulation studies and 
estimates of the accuracies of the various methods. 

Pseudo-Range 

This method of observation 1$ the one with which GPS was Initially designed to operate. It 
will be primarily used for navigation, although If used over an extended period of time would 
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provide geodetic accuracies. The observing procedure is as follows: 

{a) The receiver correlates its generated PRN codes with those of the satellite (while also 
cojnpensating for Doppler shift of the codes). The time shift (x) necessary to make this correla- 
tion results from the range to the satellite and the receiver clock synchronization error (at). 
(Due to this synchronization error, this method is called fjseudo-range and not range.) 

(b) The value of x must then be corrected for first-crder ionospheric refraction, using 
either the t values obtained from both L-band frequencies or a one-frequency model. 

Therefore, subject to other minor error sources, the observation equation is 

r * CT » - ?1 + c At » C(Uj - u)* ♦ (Vj - v)* + (w^ - w)^]*+ c At 

where 

c is the velocity of light 

^s ^“s '^s ®*'‘^^*^'**®*^ satellite vector 

p (u V w)^, earth-fixed station vector 

Since c, u^, v , w^ are assumed to be known, there are four unknowns left, at, u, v, and w. An 
exact solution results if four satellites are observed simultaneously. This is the basis of the 
GPS navigation solution and of most of the GPS navigation receivers. 

Simulation studies by Fell [1980] show that using this method in the dynamic (point position- 
ing) mode, geocentric coordinate? can be provided with an accuracy of approximately 50 to 80 cm 
after one day of observation, or 25 to 40 cm after five days of observation (see Table 1). Note 

Table 1 Effect of Systematic and Random Error Sources on Dynamic 
Point Positioning Using One-Hour Satellite Tracking 
Intervals 
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from [Fell, 1980, Table 6.2.1] 

that this type of observation does not require observations to four satellites simultaneously, 
thereby simplifying the receiver and allowing continuo' s use anywhere in the world (since one sat- 
ellite will always be well above the horizon). However, better oscillators (clocks) are assumed 
to have been used than those of the navigation receivers, and tropospheric refraction must be more 
carefully modeled. If two receivers are used to simultaneously observe the same satellites, high- 
er accuracy may be obtained for the' basel ine. The simulations show that accuracies on the order 
of 3 cm for baseline lengths are possible, with tropospheric refraction and the receiver clocks 
introducing the most error (see Table 2). 
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Table 2 Effect of Systematic and Random Error Sources on Baseline 
Determination Using One<>Hour Satellite Tracking Intervals 
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from [Fell. 1980, Table 6.3.1] 
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Doppler Observations 


The observation method for Doppler Is very similar to that of the pseudo* range. As explain- 
ed for pseudo-ranges, a receiver must receive and correlate the PRN and 0 codes. PoUoMlng this, 
the codes are removed, the resulting continuous wave carrier Is compared with the local (receiv- 
er) oscillator, and the beat frequency Is measured to obtain the Doppler count. This Is the 
method presently used In the TRANSIT satellite system. Alternatively, the receiver must continu- 
ously shift Its generated PRN codes to match the satellite's codes to correct for the Doppler 
effect. A measure of this shift also corresponds to the Doppler count With either method, 

the counts are then corrected and used In the standard range difference Doppler observation 
equation: 

^jk • 


where 



Is the reconstructed satellite frequency 
is the local generated frequency 

are the end and start times of the count at station J 


(The above equation neglects systematic error effects such as oscillator drift or tropospheric 
refraction.) According to Table 1, the accuracy of geocentric positions obtained by this method 
appears to be between 65 to 150 cm after one day of observation or 30 to 70 cm after five days of 
observation. This again assumes tracking of Individual satellites for one hour at a time, with 
most of the uncertainty due to the satellite ephemerls. If a baseline Is measured (see Table 2), 
accuracies of about 5 cm are also obtainable, assuming that receiver white noise contributes 
about 1 cm error. (This may be too low. The true value will depend on the actual Instrument In 
use.) The limiting factor seems to be the receiver white noise. 


Interferometric Observations 

Basically three different types of Interferometric observations have been proposed. Before 
explaining them individually. It Is appropriate to note what they have In common. All of the 
Interferometric methods would be used for baseline (actually baseline component) determination 
only. Thev assume that two stations will be observing the same satellite pass. This Implies 
that the satelllte(s) must be visible to both stations, so that for the GPS satellites, there Is 
a limit of about 4000 km on the length of the baseline to be measured, and that the data at both 
stations must be correlaced for a solution. One method makes use of measurements of the phases 
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of what are called the "reconstructed carrier" si 9 na 1 s« with knowledge of the original modula- 
tion code being used in generating these signals. A second method makes use of knowledge of the 
general structure of the code, but not the code itself. A third method is similar to astronomi- 
cal radio interferometry, but uses the GPS signals as the "noise" sources. No knowledge of the 
signal structure or the code is needed with this method. These methods are discussed below [NRC, 
1981]. 


The iiMt method makes use of nearly simultaneous measurements to two or more satellites 
from each ground station; the effects of station frequency standard instabilities can be almost 
eliminated in any of the methods. In the most favorable cases, the main accuracy limitation ap- 
pears likely to result from the uncertainties in the radio wave propagation corrections due to 
water vapor content of the atmosphere. Water vapor radiometers probably will need to be used at 
both ends of the baselines to infer the integrated water vapor content along the lines of sight 
to the satellites to achieve centimeter-level accuracy, unless the baseline lengths are short. 

This is because the atmospheric water vapor content is likely to be quite variable in time and 
somewhat inhomogeneous spatially at many locations. 

The basic principle of using the measured difference in phase for signals received at two 
sites to determine the component of the baseline in the direction of the source is well known 
from astronomical radio interferometry. Such measurements would be simplified if the signals 
transmitted by the GPS satellite were sinusoidal. Unfortunately, instead of being nearly mono- 
chromatic, the GPS signals have strongly suppressed carriers. The code modulation produces a 
spectrum spread over about 20 Mhz for each of the two frequency bands transmitted. However, sig- 
nals equivalent to the carriers can be reconstructed in the ground receiver without much loss in 
the signal -to-noise ratios if the code used in the satellite is known. The basic idea is to gen- 
erate a local oscillator signal that is 180" phase modulated in the same way as the signal from 
the satellite of interest, so that the beat betwee" the two signals is nearly free of the modula- 
tion effects. The beat is used to produce a clean phase-locked output signal which will change 
its phase by one cycle each time the radio path length to the satellite changes by one wavelength. 
This is called a reconstructed carrier signal. 

The reconstructed carrier phase method [Bossier et a1., 1980] makes use of nearly simultane- 
ous phase measurements for a particular satellite by receivers at both ends of a baseline. Only 
the phase difference between the signals received at the two ends is needed in the analysis, so 
that satellite clock instabilities cancel out. If the two receivers observe the phase differ- 
ences between the signals from two satellites simultaneously (double-phase differencing) both 
satellite and ground dock errors cancel out. For example both receivers might switch simultane- 
ously every second to a different satellite. In this case short switching time is needed to 
reduce the stability requirements on the receiver clocks. 

Fell [1980] has performed simulation studies using the interferometric double phase differ- 
encing method also. Assuning observations on pairs of satellites from two stations over six 
hours, he determined that baseline errors were on the order of 2 to 4 cm for 100 km baselines 
(see Table 2). Other simulations show that the observation time may be shortened considerably 
[Bossier et al., 1980]. 

The second method is similar to the one just discussed in that signals equivalent to thn re- 
constructed carrier signals are generated in the receivers. However, this is done without kxact 
knowledge of the code being required. 
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Tmo versions of this method are being pursued currently. At the Jet Propulsion Laboratory, 
the SERIES system which Is being built aod tested uses directional antennas about 1.5 m In diam- 
eter to track Individual GPS satellites (Fig. 1). Knowledge of the approximate switching rate 
for the code is employed to resolve ambiguities without the requirement for a long observing 



Fig. 1 The SERIES Directional Antenna and Van 
(Photo coufUtiy oi JPL,) 

time at each site, but no other Information on the code Is needed [MacDoran, 1981]. Another ver- 
sion of this method [Counselman, 1981] uses small, nearly omnidirectional antennas (Fig. 2) and a 
different data processing technique. Results so far Indicate that reflection of the satellite sig 
nals from the ground or nearby structures Is not a serious problem, even with the nearly omnidi- 
rectional antennas. 



Fig. 2 The MITES Nearly Omnidirectional Antenna Placed on 
tne Ground (Photo oi C.C. Counae/mon III) 
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Th< ttuAd met/iod Is analogous to the one used with radio astronomical signals. The use of 
1ong»basel1ne radio Interferometry with noise transmitted from a satellite was actually dam* 
onstrated some time ago. An adaptation of this method for accurate baseline determinations with 
the spread* spectrum GPS signals as the radio noise sources was the original SCRIES proposed by 
MacOoran [1979, 1980]. Since the received flux density Is about 10* times larger than for extra* 
galactic radio sources, the apparatus can be much less expensive and more compact than that needed 
for astronomical long*basel1ne Interferometry. 

The second method discussed above appears to be superior to the third one because It requires 
less data recording and processing to achieve a given s1gna1*to*no1se ratio. However, a comparl* 
son of the advantages of the first and second methods Is difficult because of uncertainty concern* 
ing the future availability of the code to non*Department of Defense users. The Department of 
Defense Is reported to be considering adopting a policy under which Incomplete knowledge of the 
characteristics of the GPS signals by nonapproved users would limit the accuracy of absolute posi* 
tion determinations, but not of relative position determinations for which real time results are 
unnecessary. However, how this policy would be iinplemented Is not known at present. It would be 
highly valuable If Information on the future availability of the code and other related questions 
could be provided soon. 

An additional Interferometric approach which has been investigated jointly by MIT and Draper 
Laboratories Involves placing some small supplemental transmitters on future GPS satellites 
[Counselman and Shapiro, 1979]. This would permit con-lderable simplification in GPS receivers 
for geodetic uses, which is highly desirable and would avoid ambiguity problems. However. It is 
not known what the chances are of adding the nwcessary equipment to the satellites. In the 
interim, demonstrations of both the reconstructed carrier interferometry method and the spread* 
spectrum Interferometry method are being carried out to assess the accuracy likely to be achieved 
with the proposed future system. 

Although there ar> still many uncertainties, the prospects seem good that all of the GPS In- 
terferometric methods being developed will achieve 3 cm accuracy for the three baseline components 
at a high confidence level. For very short baselines, about 100 m in length, Counselman [1981] has 
demonstrated 1 cm accuracy In each of the components. However, a much better knowledge of the 
range of errors associateci with tropospheric water vapor In different regions Is needed In order to 
clarify the prospects further. It seems clear that at least two of the GPS interferometric methods 
will be preferred over the u«e of astronomical long*basel1ne Interferometry for baselines which are 
short enough so that the error contributions due to the GPS orbit uncertainties do not become sub- 
stantial. Low GPS satellite orbit uncertainties are expected when accurate interferometric track- 
ing data is available from a well -distributed set of fixed ground stations [MacOoran, 1379; Coun- 
selman and Shapiro, 1979]. 

Summary of Instruments 

Table 3 lists geodetic GPS receivers, their stage of development, and different characteris- 
tics. This list is not meant to be complete, but merely gives some Indication of the types of re- 
ceivers presently being worked on. 

It is obvious that with the current rate of development, geodetic receivers will be readily 
available by the time the GPS system Is in full operation, with the interferometric instruments by 
then providing the highest b.itse11ne accuracies. However, the pseudo-range and Doppler Instruments 
are already available and capable of providing station positions with accuracies similar to those 
obtainable using the TRANSIT system of satellites In conJunct1o.h with the by now "classical" Dop- 
pler integrator Instruments. 
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Table 3 GPS Receivers (User Sepentt Geodetic Use Only) 


Name 

CNA-782 

Oroaniaatlon 

Canadian 
Marconi Co. 

latki 

Available 

commercially 

ill! 

Semi -portable 
(10" X 17" 
antenna) 

Observable 

Pseudo-range 
Doppler phase 

Reference 
Blaha 1980 

ST1-5010 

Stanford 

Telecommunications 

Inc. 

Available 
coneercial ly 

Van mounted 

Pseudo-range 

Doppler 

Perreault 1980 


Hagnavox 

Under devel- 
opment 

Portable 


Deem 1980 


Texas 

Instruments 

Under devel- 
opment 



Ward 1980 

SERIES 

Jet Propulsion 
Laboratory 

Testing 

Van mounted 
(1.5 m dish 
antenna) 

Phase differ- 
encing 
Range 

MacOoran 1981 

MITES 

Massachusetts 
Institute of 
Technology 

Under devel- 
opment. antenna 
being tested 

Portable 

Phase differ- 
encing 

Counselman 1981 
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September 21, 1981 


From: Ivan I. Mueller 

Subject: Minutes of CSTG Meetings, Munich, September 2, 1981 


The meeting of the Commission on International Coordination of Space 
Techniques for Geodesy and Geodynamics (CSTG) was opened at 8:35 AM by the 
President, Ivan I. Mueller. Kenneth I. Daugherty was asked to act as secre- 
tary. The agenda for the meeting was distributed (Enclosure 1), items 1 through 
7 being covered in this meeting. Sixteen people were present either as observ- 
ers or members of CSTG (Enel. 2). 

1. Pres1dent*s Report . Presented by Mueller. In response to questions 
regarding the MEDOC project, C. Boucher stated that data collection has been 
suspended on MEDOC until some Improved methods for data analysis are available. 

2. Subcommission Report . Boucher presented the report on the Sub- 
commission on Standards. Tnis report, together with the reports of the other 
subconimissions, is found In CSTG Bulletin No. 3; 

3. Participation In Tokyo . Mueller reported that a half day (AM, May 14, 

1982) had been set aside for CSTG at the LAG meeting in Tokyo In 1982. Present 

plan calls only for a business session, but there is some pressure to present 
some scientific papers. Robbins suggested a report on Project MERIT. Ch. Reig- 
ber suggested review papers on projects related to CSTG to cover scientific back- 
ground and status. Mueller asked who would be in Tokyo, and the following re- 
sponded: Robbins, Boucher, Giacaglla, Relgber and Prilepin. It was agreed that 

review papers of the above type would be solicited by the President. 

4. Short Course in Manila . The Philippines proposed hosting a one-week 
short course on the applications of space techniques for geodesy and geodynamics 
to be held in Manila. It Is proposed that If a special program were to be 
established, it should follow the I AG meeting in Tokyo In 1982, with from 10 to 
12 experts going to Manila to offer technical lectures. CSTG, through COSPAR, 
has requested funding from ICSU for the travel. The President requested that 

ACTION everyone think about this and participate if possible. Let Mueller know in the 
next two months. 
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5. Project ALGEDOP . Mueller reviewed a letter from A. Marussi propos- 
ing ALGEDOP (Alpine fieoid Doppler Project) (Enel. 3). The participants expressed 
support for this project provided that it can be soundly justified on a scientific 
basis. In this regard the main questions were the need for such a geoid and 
whether the accuracies achievable through Doppler are really superior to those 
obtainable through gravimetry. 

6. WARC/SFCG Frequency Allocation . Reigber presented details of POPSAT. 

To support this project^thi^e is a need for an up/ down link frequency allocation. 
The desired frequency is in the 2 gHz range. The POPSAT group would like CSTG 
to appeal to the appropriate authorities in support of setting aside the re- 
Quired frequency for their system. At this point, action on the item was suspended 
and was resumed during the afternoon session at 4:35. At that time Reigber briefly 
reviewed the international frequency allocation structure and procedures. A reso- 
lution was passed that the President write to the executive secretary of the 
International Space Frequency Coordination Group (SFCG) putting forth support for 
the frequency allocation for POPSAT (Enel. 4). 

7. Other Business . As there was no other business the session was 
adjourned to be reconvened at 2:30 for a joint session with lAG Commissions X 
and XI. 


Joint meeting of CSTG with lAG Commissions X and XI convened at 2:30, 
September 2, 1981. Mueller, R.O. Coker for Commission XI and R. Sigl for 
Commission X jointly chaired the meeting. See Enel. 1 for the agenda. Items 
8 through 16 were covered during the meeting. Twenty-six persons attended 
(Enel. 5). 

8. Introduction . Coker opened the meeting by reviewing the background 
of the Doppler project for Africa. He requested Mueller to continue the meeting 
and to cover the items on the agenda. 

9. Resolutions of the Addis Ababa Symposium . Mueller gave the history 
of CSTG activities related to the Doppler project for Africa and outlined the 
sister nation concept which is being advocated by CSTG and Commission XI (Enel. 


10. Reports on Existing Projects or Proposals . Mueller called on the 
various representatives present to report on existing projects or proposals. 
The first of these was the Upper Volta project. P. Richardus reported on the 
Dutch project. It consists of 16 Doppler points as part of an overall survey 
program. The observations are being done with Marconi equipment. 

L. Sjdberg reported on the Tanzania and Zambia projects proposed by 
Sweden (Enel. 7). 

Daugherty reported on DMA projects and proposals. DMA has done work in 
Africa for a number of years and has established 60-80 points by Doppler point 
positioning. DMA currently has agreements with Egypt and Sudan which will 
involve additional Doppler points. DMA would like to establish bilateral 
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agreements with other African nations. DMA would furnish Doppler receivers 
and full teams (two persons) and would pay all expenses of their teams. Sister 
countries would be expected to supply local support and participate In the field 
program. All data would be reduced by DMA and results furnished to the sister 
nation for release. DMA 1s also proposing to serve as a computation center for 
the African Doppler project to reduce data using the precise ephemerls to estab- 
lish Doppler point positions. 

E. Reinhart reported on the Ivory Coast project (Enel 8). 

S. Krynski described related activities by Poland (Enel. 9). Based on the 
Polish report Mueller asked the Polish representatives If they would prepare a 
plan for the Doppler survey of Africa to be presented at the Nairobi meeting In 

ACTION November and send It to CSTG and to Commissions X and XI by Oc*tober 15. The 
Polish representatives agreed to do so. 

11. Results of the Surveys by Questionnaires . Mueller reviewed the results 
of the questionnaire which he had sent to countries who might be able to give sup- 
port for the Doppler survey of Africa. A summary based on these questionnaires Is 
provided In Enel. 10. This summary Indicates that there Is Interest in at least 

14 countries which together own more than 110 Doppler receivers of various types. 
Subsequent to the meeting Bulgaria and Brazil also forwarded their questionnaires 
which indicate interest in participation by a total of three Doppler receivers 
(one MX, JMR-1, CMA 761). A proposal for selection from these will be presented 
at the Nairobi meeting. 

12. PeyelopTOnt of Standards . Mueller introduced copies of the new 
Canadian manual on Doppler surveys and two DMA manuals on the same subject. He 
proposed that a committee of experts chaired by Boucher review the documents and 

ACTION make a recommendation by October 15 on standards of survey and documentation to be 
used In the African campaign, Boucher selected B.A. Sikllo, W. Schlueter and 0. 
Critchley as members of his subgroup. The proposal for standards will be presented 
In Nairobi. 

13. Regional Centre for Services In Surveying and Mapping In Nairobi . 
Mueller stated that CSTG would encourage all computations to be done at two 
computing centers. The Regional Centre for Services In Surveying and Mapping 
in Nairobi has volunteered to be the African computing center. Those countries 
which have an Interest in acting as a computing center should send a letter of 

ACTION Interest by October 25 to CSTG and Commission XI describing their capability 
and experience In Doppler geodetic computations. 

14. Where Do We Go from Here? It is important that the plans and as 
many arrangements as possible be ready before the Nairobi meeting. Interested 
parties should go to Nairobi prepared to discuss details of the African Doppler 
program. It was agreed that the President should again contact all African 
countries, inform them about the outcome of this meeting, and encourage them 

to participate at the Nairobi irteeting. 

15. The Meeting in N ovember, 1981 . At the meeting In Nairobi 

it Is the desire of the CStG and Commission XI members present that our business 
be conducted during the first week. If possible. Coker agreed to arrange for such 
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meetings on November 10 and 13. In the meantime, J.D. 0be1 , Secretary of 
the Nairobi conference Indicated that the CSTG meetings would be scheduled 
November 13, 16 or 17. CSTG Is now planning to meet on the 13th and 16th. 

As described above, during this meeting the following will be presented: 

— plan of Doppler stations 

— draft of specifications for the survey 

— designation of computing centers 

— a draft statrment which would need to be part of any bilateral 

agreement without which the survey will not qualify as 
part of this project 

— proposal for the designation of sister countries 

-- dates for the campaign (see below) 

16. Other Business . There being no other business the joint meetings 
were adjourneBT At a later Informal meeting, the suggestion was made that the 
Doppler survey should be conducted during a specified one- to three-month 
period during which time precise ephemerldes of all NNSS satellites would be 
available, sometime during 1982/83. Daugherty was requested to report at the 
Nairobi meeting on such time periods. 


Note: The enclosures, omitted here, may be obtained on request 
from the President of CSTG. 



